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("^^ Abstract. We consider the Glauber dynamics for the Ising model with "+" boundary condi- 

rvj tions, at zero temperature or at temperature which goes to zero with the system size (hence 

-. the quotation marks in the title). In dimension d = 3 we prove that an initial domain of lin- 

JjT^ ear size L of "— " spins disappears within a time t+ which is at most L^(logL)'^ and at least 

ry\ L^/(c log I/), for some c > 0. The proof of the upper bound proceeds via comparison with an 

auxiliary dynamics which mimics the motion by mean curvature that is expected to describe, 

\l on large time-scales, the evolution of the interface between "+" and "— " domains. The analysis 

' of the auxiliary dynamics requires recent results on the fluctuations of the height function asso- 

i_^ ciated to dimer coverings of the infinite honeycomb lattice. Our result, apart from the spurious 

r^ logarithmic factor, is the first rigorous confirmation of the expected behavior t+ ~ const x L^, 

f^ I conjectured on heuristic grounds [12, 6]. In dimension d = 2, t+ can be shown to be of order 

J-i L^ without logarithmic corrections: the upper bound was proven in [7] and here we provide the 

'^j lower bound. For d = 2, we also prove that the spectral gap of the generator behaves like c/L 

C^ for L large, as conjectured in [2]. 
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A long standing open problem in the mathematical analysis of the stochastic Ising model [16] 
can be described as follows. 



o 

(^ Consider the standard ±1 spin Ising model at inverse temperature /3 and zero external mag- 

^~~' netic field in a cubic box A C Z , d > 2, of side L, with homogeneous, e.g. "plus", boundary 

J> conditions outside A (see Section 2 for the precise definition). Denote by vr^ the corresponding 

Gibbs measure and assume that /3 is larger (or much larger) than the critical value /3c where the 

"plus" and the "minus" phases start to coexist. On the spin configuration space f^A consider 



^ also the continuous time Markov chain, reversible w.r.t. the Gibbs measure vrj, in which each 

spin cj^;, X G A, with rate one chooses a new value from ±1 with probabilities given by the 
conditional Gibbs measure given the current values of the spins outside the site x. Such a chain 
is known in the literature as Gibbs sampler or Glauber chain and, because of reversibility, its 
unique stationary distribution coincides with ttJ. Let fi^ denote the distribution of the chain at 
time t when the initial configuration at time t = is a G { — 1, 1} and let Tmix be the minimum 
time such that the variation distance between fi^ and vr^ is smaller than, e.g., l/(2e) for any a. 
Because A is a finite set and the chain is ergodic, the Perron-Frobenius theorem guarantees that 
Tmix is finite and it is not too difficult to prove that Tmix < e'^^^ for some constant c (see e.g. 
[16]). Despite the fact that the above upper bound is known to be the correct order of growth 
of Tjiiix when the boundary conditions are absent, i.e. free (see e.g. [24, 16]), it is known that the 
presence of the homogeneous "plus" boundary conditions drastically modifies the whole process 
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of relaxation to equilibrium. In this case the conjectured correct growth of Tmix as a function 
of L is given by the Lifshitz's law Tmix = 0{L?) (see [12] and [6]). 

Unfortunately, any polynomial upper bound on Tmix has escaped so far a rigorous analysis and 
the best known result, following a recent breakthrough [18], is confined to the two dimensional 
case d = 2 and it is of the form Tmix < exp(clog(L)^) for any (3 > /3c (see [14]). 

The conjectured L^ growth of Tmix is related to the shrinking of a "spherical" bubble of the 
"minus" phase under the influence of the "plus" boundary conditions. On a macroscopic scale 
one expects [6, 23] that the dynamics of the bubble follows a motion by mean curvature which, 
if true, implies immediately that in a time 0{Lp') the bubble disappears and equilibrium is 
achieved. 

Some of the above questions and in particular the validity of the Lifshitz's law can be formu- 
lated (and its mathematical justification remains highly non-trivial) even at zero temperature 
(/3 = +oo), a situation that was considered in great detail in [7] (see also [5], [23]). Similarly one 
can consider, as we do in the sequel, an almost zero temperature, i.e. (3 = /3{L) and increasing 
so fast with L that thermal fluctuations become irrelevant on the relaxation process. In the 
extreme case (3 = +oo, the stationary distribution is concentrated on the "plus" configuration 
(i.e. all spins equal to +1) and the Glauber chain evolves towards it without ever increasing the 
spin energy. The mixing time Tmix becomes closely related to the hitting time of the "plus" 
configuration starting from all minuses and it is not too hard to prove, by induction on the 
dimension d (starting from the case d = 2, which requires non-trivial work [7]), an upper bound 
of the form Tmix ^ cL'^. However the inductive argument, which in d = 3 simply boils down 
to comparing the true evolution of the original cubic bubble of minuses with an auxiliary chain 
in which the L two-dimensional square layers of the bubble disappear one after the other start- 
ing e.g. from the top one, completely neglects the interesting cooperative effect in which the 
whole bubble, starting from the corners, is eroded by a "mean curvature effect" which enhances 
considerably its shrinking. In other words, the Lifshitz's law in d = 3 cannot be proved or 
even approached closely without considering in detail how a two-dimensional curved interface 
separating the pluses from the minuses evolves in time. 

The first main contribution of the present work is a strategy to attack and solve the above 
problem at zero temperature for d = 3 (with logarithmic corrections) and for d = 2 (with different 
constants in the upper and lower bounds). We refer to Section 3 for the precise statements. In 
the challenging three dimensional case our method involves the use and adaptation to our specific 
situation of two different and beautiful sets of results concerning: 

(i) the Gaussian Free Field-like equilibrium fluctuations of random monotone surfaces and 
their connection with random dinier coverings of the two-dimensional hexagonal lattice 
(see [4, 8, 9, 11, 21]); 
(ii) the mixing time of a Glauber chain for monotone surfaces [25] . 

The key point of our approach is to prove that the evolution of the Glauber chain is dominated 
by that of another effective chain which follows a motion by mean curvature, only slowed down 
by a logarithmic factor (in the radius of the bubble) . In turn the evolution of the bounding chain 
is constructed by "peeling off" at time t a layer of logarithmic width from the bubble of radius 
Rf. The peeling process is realized by letting relax to equilibrium a mesoscopic spherical cap of 
radius 0{y/Rt x polylog{Rt)) and height O {poly log (Rt)) (where polylog{x) stands for a suitable 
polynomial of log x) centered at each point of the surface of the bubble. Results (i) above are 
essential in order to prove the peeling effect while the results (ii) prove that the overall effect 
occurs on the correct time scale 0{Rt x polylog{Rt)) ■ We strongly believe that our approach 
can be helpful in solving other related problems. 



The second contribution, this time for the model in dimension d = 2, are upper and lower 
bounds, linear in L~^ and uniform in /3 > clogL, c large enough, on the spectral gap (see (2.6)) 
of the Glauber chain. In [2] for d = 2 the upper bound was proved (apart from logarithmic 
corrections and with constants depending on (3) for any 13 > Pc and it was conjectured to be the 
correct behavior of the spectral gap (to be compared with the L^ scaling of the mixing time) . In 
our case the proof of the lower bound (the most interesting one) has an analytic flavor. We first 
unitarily transform the original Markov generator of the Glauber chain acting on i^{Q.j\^,T:^) 
into a new matrix acting on £^(i7A) with the flat (i.e. counting) measure, whose off-diagonal 
elements corresponding to spin transitions which do not conserve the energy vanish very fast as 
j3 —7- oo. Such a property, for /3 > clogL, allows us to write the new matrix into a block-form 
plus a remainder whose norm is very small with L. Since each block describes a suitably killed 
Glauber chain (killing occurs as soon as the spin energy decreases) the desired bound follows by 
a probabilistic analysis of the killing time for each block. 
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2. Model and preliminaries 

Given x,y G Z , we write x ~ y if x is a nearest neighbor of y. To each point x G Z is 
assigned a spin ax which takes values in { — 1, -|-1}. Given A C Z , we let 

(9A := {x G Z'^ \ A : 3y G A such that x ~ y} 
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and r^A := {— 1,+1} be the set of all possible configurations in A. The Gibbs measure in a 
finite domain A for a boundary condition rj £ ^dA at inverse temperature /3 > is given by 

<i^) = -^^ (2-1) 



with 



^a(^) •= - X] ^^^J' ~ Z] ^^^J' (^-^^ 



x,y<^K x€A,y€dA 

x^y x^y 

and 

Zl := ^ e-f'^Ai-l (2.3) 

When no danger of confusion arises, we will write just vr for vr^. We will consider also the case 
P = +00, in which case ir^ is taken to be the uniform measure over all configurations a G Q\ 
which minimize H^j!^. 

The dynamics we consider is the Glauber dynamics {o'^{t)}t>o, where £, = cj^(O) is the initial 
condition. To each x G A is associated an independent Poisson clock of rate 1. When the 
clock labeled x rings, say at time s, one replaces ax with a value sampled from the probability 
distribution vr^^o-«{s)(")) where 

-n-xA-) ■='^(.-Wy^y 7^^)- (2-4) 

We denote the law of cr^{t) for a given time t as /U^. It is well known that {cr^(t)}i>o is a Markov 
process, reversible with respect to the equilibrium measure vr. 

We are interested in the order of magnitude of the "mixing time", defined for some e G (0, 1) 
as 

Tmix(e) = inf{t > : sup ||^^ - 7r|| < e} (2.5) 

(tGSIa 



where 



^11 = 9 IZ l^(^)-'^(^)l 






denotes the total variation distance between two probability measures. When e = l/(2e), we 
will just write Tmix for Tmix(e). Another quantity we will focus on is the spectral gap, 

, . ^ <(/(-£)/) 
gap = gapX = mf — — — 2.6 

Var^^(/) 

where Var^»7 denotes the variance w.r.t. tt^, £ is the generator of the dynamics, 

Cf{a) = Y,{7Tx,M)-ficr)) (2.7) 

xeA 

(a more explicit expression for the generator is given in Section 9.1) and the infimum is taken 
over non-constant functions / : Qa i— )■ M. The inverse of the gap (respectively the mixing time) 
measures the speed of convergence to equilibrium in L'^{diT) (resp. in total variation norm). 
Also, it is well known that 

sup ll/u^ -vr|| < e^Li/TmixJ^ (2.8) 

(tGSIa 

which actually holds for any reversible Markov process. 



2.1. Preliminary results. 

Monotonicity and global coupling. For notational clarity, in this subsection we indicate explicitly 
the dependence on the boundary condition rj in /i^ , the law of the dynamics <T^{t) at time t, ^ 
being the initial condition. The Glauber dynamics enjoys the following well known monotonicity 
property. One can introduce a partial order in Q\ by saying that a < a' if ax < o'^. for every 
X G A. Then, one has 

4,, ^4Vif e<r,r?<V (2.9) 

where -< denotes stochastic domination (one writes ^x < v \i fJ,{f) < ^(/) for every increas- 
ing function /, i.e. f{cr) < f{cr') whenever a < a'; an event will be called increasing if its 
characteristic function is increasing). In particular, letting t — t- oo one obtains 

ttI ^ nf. (2.10) 

Also, it is possible to realize on the same probability space the trajectories of the Markov chain 
corresponding to distinct initial conditions ^ and/or distinct boundary conditions ij, in such a 
way that, with probability one, 

<yf,{t) < (J^^,{t) for every t > 0, if ^ < ^',r/ < r?'. (2.11) 

This will be referred to as the "monotone global coupling" : its law will be denoted by P and the 
corresponding expectation by E. 

Throughout the paper we will apply several times the above monotonicity properties: for 
brevity, we will simply say "by monotonicity..." instead of referring explicitly to (2.9)-(2.11). 

Comparing /3 = +oo and (3 > ClogL. In this section, it is understood that A is the cubic 
domain 

kL:={-L,...,LY (2.12) 

of side 2L + 1 and that rjx = + for every x G ^A^, so that we omit A and r/ from the notations. 
In this work, we consider the situation where /3 grows with L, and in particular 

/3g (ClogL, oo] (2.13) 

for a sufficiently large C. 

Note that, for /3 = +oo, the equilibrium measure concentrates on the all "+" configuration: 
TT°^{a) = lo-=+- Given c > 0, one can choose C large enough so that 

lk-7r°^||<— (2.14) 

for every (3 in the range (2.13) (this follows from easy Peierls estimates). Moreover, for every 
initial condition o" G fi, one can find a coupling between the dynamics at /3 = +oo and /3 > 
ClogL such that they coincide until time L'^ with probability at least 1 — 1/L^ (just compare 
the transition rates). In particular, this implies that 

\K - ih''^\\ < j-c ^°^ ^^^^y * - ^^- (2.15) 

We will see later (cf. Remark 1) that, as an immediate consequence of Theorem 1, one actually 
has the estimate 

lim sup sup ||;U^ — ;U^'°^|| = 0. (2.16) 
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We define the random time r+ as the first time when all spins are "+" , starting from the all 
"— " configuration: 

T+ := inf{t > : cr^(t) = +1 for every x G Al}. (2.17) 

Observe that, under the global coupling, if /3 = +oo, then ai{t) = +1 for every t > r+ and for 
every initial condition ^. 

Notational conventions 

• For notational clarity, quantities referring to the (3 = +oo dynamics will have a super- 
script oo (we will write for instance //f' °° and T^;^), while (for lightness of notation) we 
will not put a superscript /3 when /3 < oo; 

• our main focus will be on the case of the all "+" boundary condition {rjx = +1 for every 
X G dh) and we will just write 77 = + in this case; 

• we let P denote the law of the process {o"~(t)}f>o at /3 = +00, started from the "— " 
configuration, with boundary condition r/ = +; 

• given a G f^A and x G A, we will denote by o"^^^ G J^a the configuration obtained by 
flipping the spin at x to —ax- 

Based on the observations (2.14) and (2.15), one can use the time r+ to estimate the mixing 
time for (3 large: 

Lemma 1. For every e G (0, 1) one has 

T-,(6) = inf{t > : P(t+ > t) < e}, (2.18) 

where we recall that P is the law of the process {(y{t)}t>Q for f3 = +00, with "+" boundary 
conditions. 

Moreover, fix e & (0, 1) and 6 > and assume that, for some cq > 0, T^j^(e) < L*^" for L 
sufficiently large. If (3 £ (ClogL, +00) with C sufficiently large, then for every L large one has 

T„,ix(e + <5)<T-,(e) (2.19) 

and 

T^ix(e-<J)>T~,(e) (2.20) 

Proof of Lemma 1 . Recah that 7r°^(o-) = 1^=+ so that ||/if'°^ - 7r°^|| = 1 - fif'^ia = +) < 
1 — P(cr^(t) = +) (where the inequality follows from monotonicity) and (2.18) is immediate. 

To prove (2.19)-(2.20), it is sufficient to choose C sufficiently large so that (2.14) and (2.15) 
hold for c = Co and to take L sufficiently large so that 2/L'^° < a. 

Lemma 1 

D 

3. Results 

3.1. Three-dimensional model. Consider the 3D Ising model in the cubic box A^ = {—L, . . . ,L} 
with boundary condition rj = + and (3 = +00. The main result of this work is: 

Theorem 1. There exists a positive constant c such that for L > 2 one has 



3 



clogL 



<T+ < cL^ (log L)^ 



>l-f (3.1) 



Remark 1. 



(1) As we mentioned in the introduction, the previously known upper bound for r^. was of 
order L^ [7, Th. 1.3], while (3.1) matches the heuristically expected behavior, except for 
the logarithmic factor. On the other hand, we are not aware of previously known lower 
bounds, except for the trivial estimate r+ > cL; 

(2) via Lemma 1, once we prove (3.1) we also have that 

n 7 < Tmix < cL\logLy (3.2) 

for f3 > ClogL and C large; 

(3) To prove (2.16), observe first of all that 

II O" CTiOOll ^ II cr II ( II cr.oo OOII 1 II OOII /o o\ 

Wfit - IJ't II < ll/^t -^11 + ll/^t -^ ll + lk-vr ||. (3.3) 

The first two terms are estimated through (2.8), taking say t>L^ and using the upper 
bound (3.2) for the mixing time. The proof is concluded thanks to (2.14)-(2.15) (if C is 
chosen such that these bounds hold for c = 3). 

3.2. Two-dimensional model. For the 2D Ising model in the cubic box Al = {~L., . . . , L]^ 
with boundary condition r/ = + and /3 = +oo we have 

Theorem 2. There exist positive constants c, 7 such that 

P (cL2 < r+ < (1/c)l2) > 1 _ e'^^. (3.4) 

for every L > 1. 

Remark 2. Thanks to Lemma 1, this implies that ciL^ < Tmix ^ C2L^ for /3 > ClogL. We 
mention that the bound P(t+ > (l/c)L^) < exp(— 7L) was proven in [7, Th. 1.3]; as for the 
lower bound for r+, the authors of [7] proved the following weaker result: for every 6 > 0, 
P(r+ < L^/(logL)^''' ) tends to zero as L — )• 00. For a simplified version of the dynamics, the 
authors of [5] proved the sharp behavior t+ '^ const x L^. 

In two dimensions, always with "+" boundary conditions, we also have sharp bounds on the 
spectral gap: 

Theorem 3. There exist positive constants C, c such that for every (3 G [ClogL, +00), one has 

|<gap<^. (3.5) 

The reason why we require /3 < 00 is that for /3 = +00 the equilibrium measure is concentrated 
on a single configuration (the all "+" configuration), and the spectral gap (2.6) is not well defined 
(all functions have zero variance with respect to vr). 

Remark 3. On the basis of the analysis of a one-dimensional toy model of evolution of a bubble 
of "— " phase inside the "+" phase, it was conjectured in [2, Sec. 7] that in dimension d = 2 the 
spectral gap behaves like const/ L for every (3 > f3c- Our result (3.5) is the first confirmation 
of this conjecture, in the /3 — t- 00 limit. Let us remark also that the authors of [2] showed, via 
the construction of a suitable test function, that gap < ci{l3){\ogLY^^"' / L for every /? > /3c. 
However, when /3 grows with L as in (2.13), the constants ci^2(/3) possibly diverge. Therefore, 
even the upper bound in (3.5) is a qualitative improvement over known results. 

The proof that gap > c/L for (3 > C log L extends easily to the three-dimensional model (cf. 
Section 9.1) but, as we discuss in a moment, there is no reason to believe that this is the correct 
behavior of the spectral gap in d = 3. 

Two interesting questions which are left open by the above results are the following: 
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(1) What is the order of magnitude of the spectral gap in dimension d = 3? One expects it 
to be much larger than in dimension d = 2, possibly of order (logL)"'^ for some c > 0, 
or even of order 1 [2, Sec. 7], [6]; 

(2) How to remove the logarithmic corrections in the upper and lower bounds for Tmix in 
dimension d = 37 A suggestive indication that indeed Tmix > cL^ comes from the 
following argument. If one defines the "entropy constant" (or modified Log-Sobolev 
constant) as 

vr(/log/) 

Cent ■■= sup —- (3.6) 

7r(log /(-£)/) 

where the supremum is taken over positive functions / which verify TT{f) = 1, it is 
possible to exhibit a test function which gives Cent ^ cL^ for some constant c which 
is positive, uniformly for /3 > ClogL with C large. On the other hand, it has been 
conjectured [19] (and verified in various explicitly solvable examples) that for every 
reversible Markov process the mixing time is lower bounded by Cent times some universal 
constant. 

The organization of the paper is as follows. Theorem 1 is proven in Section 4 (upper bound 
on r+) and in Section 7 (lower bound). The upper bound requires precise estimates on the 
fluctuations of dimer coverings of the hexagonal (or honeycomb) lattice, see Section 5, and on 
the mixing time of a dynamics on plane partitions or monotone sets, proven in Section 6. As 
for the 2D model. Theorem 2 is proven in Section 8 and Theorem 3 in Section 9. 

More notational conventions 

• In the proof of the results, c, c', c" etc. denote positive and finite constants (independent 
of L and /3) which are not necessarily the same at each occurrence. 

• If ti G M'^, it will be understood that u^""' , a = 1, . . . ,d are its Cartesian coordinates. 

• d{-,-) will denote the Euclidean distance in M.'^. 

• dB will denote the geometric boundary of a set B CM.'^, except when B C Ij'^, in which 
case dB := {x ^ 'Z'^ \ B : 3y ^ B such that d{x, y) = 1}. 

• N denotes the set of non-zero integers {1,2,...}. 

4. The mixing time in dimension d = 3 (upper bound) 
4.1. A basic tool: dynamics of discrete monotone sets. 

Definition 1. We say that F C N^ zs a positive monotone set if z £ V implies y € V whenever 
y G N^ and y'"' < z^^',a = 1,2,3. The collection of all positive monotone sets will be denoted 
by E+, which is partially ordered with respect to inclusion. When V £ S"*" is a finite set, it is 
usually called a plane partition, the two-dimensional generalization of an ordinary partition (or 
Young diagram). 

Given V^,Vo,V^ in S+ such that V^ ^ Vq ^ T/+, we define a dynamics {V^ "}i>o on S"*" 
such that V C y/" C V+ for all times, with initial condition V^% = Vq. Let A := V+ \ V~ 
and associate to each z £ A an independent Poisson clock of rate one. When the clock labeled 
z rings at some time s: 

• if both sets Vg^z- '■= Vs \ {z} and Vs^z,+ '■= Vs U {z} belong to S^, then we replace Vg 
with Vs^z,a, with a chosen between "+" and "— " with equal probabilities 1/2; 

• otherwise, we keep Vg unchanged. 



The link with the /3 = +cx) dynamics of the Ising model is straightforward: Consider the Ising 
dynamics in the domain A := V^ \ V~ , with initial condition 



+ 


if 
if 


z eVoDA 
zeA\Vo 


if 
if 


z G (Z3 \ N3) U V 



^^ = \ ^ :, ..a\ t/. (4-1) 

and with boundary conditions 



Vz - { ;f ^ r- I'T/'i \ RJ3A , I T/- (4-2) 

and let Mt := {z G N^ : oi(t) = -}. Then, Mt has the same law as V^° . 

To avoid any risk of confusion we choose different notations for the Ising dynamics and the 
monotone set dynamics: we call v^ ° the law of V^ ° and p := py± its (reversible) invariant 
measure. Of course, py± is just the uniform measure over the positive monotone sets V G S"*" 
such that V^ ^ y C y+. The monotonicity properties of the Glauber dynamics discussed in 
Section 2.1 immediately give: 

Proposition 1. IfV^,W^ G S+ and V^ C W^ , then py± ^ pw±. 

The next key result quantifies the mixing time of the monotone set dynamics as a function of 
the shape of A: 

Theorem 4. Let 

Z):=max{d((z«,z(2)),(yW,y(2)))} 

and 

H := max{\z^^^ -y^^^l : z,y £ A and y^"^ = z^''\a = 1,2}. 
If H < D, the mixing time of the monotone set dynamics is 0{H'^ D^ (log D)^). 

We believe that the correct behavior is O(D^log-D). The spurious factor H^ is potentially 
dangerous, so we will take care to apply Theorem 4 only in situations where H is small, say of 
order of a power of logD: this is a crucial step for the proof of the upper bound (3.2), which 
differs from the expected behavior 0{Lp') only by logarithmic corrections. We give the proof of 
Theorem 4 in Section 6. 

4.2. Upper bound on the mixing time. Let 

Sr '-^ '^ n Br, 

where 

Br :={xGM^ : d(x,0) <r] 

is the ball of radius r centered at the origin. By monotonicity, the claim (3.1) follows if we prove 
the following result for the dynamics in S^l with "+" boundary conditions, started from "— ": 

Proposition 2. Consider the dynamics in S^l with "+" boundary conditions on dS^L. There 
exists c > such that for every L > 2 and t = cL^(logL)'^ one has 

P {3x G S5L \ Sl such that a^{t) = -) < c/L. (4.3) 

The reason why this result implies (3.1) is that the set 5*51, \ Sl contains a cube, call it Q, 
which is a translate of A/, and, by monotonicity, the marginal in Q of the evolution in S^l is 
stochastically dominated by the evolution in Q started from "— ", with "+" boundary conditions 
on dQ. 

In order to get Proposition 2, we prove 
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Proposition 3. Consider the dynamics in Sl with "+" boundary conditions and let x £ Sl\ 
Sl-i. There exists c > such that for every L > 2 

P [a- (t) = + for all times t G [c L(log L)^ L^]) > 1 - c/L^. (4.4) 

Proposition 2 then easily follows: 

Proof of Proposition 2 (assuming Proposition 3) By the union bound, one then deduces that 

P (3a; £ Sl\ Sl~i and t G [cL(logL)^ L^] such that o-^(t) = -) < c/L^. (4.5) 

For the dynamics in S^l we prove, by induction on i = 1, . . . , 4L, that 

P{Ai) := P (3x G S^L \ Sr,L~i and t G [5ciL(log(5L))^ L^] such that cr~{t) = -) < -^, (4.6) 

which for i = 4L implies (4.3). For i = 1, this follows from Eq. (4.5). If the claim holds for 
some 1 < i < 4:L, then 

P(A+i) < P(A) + P(A+i| A') < ^ + P(A+i| A') (4.7) 

where for any event A we let A'^ denote its complement. Next, by monotonicity, 

P(A+i| A) < P (3x G S5L-i \ S5L-i-i,t G [5c [i + l)L{log{5L)y, L^] : a,{t) = -) (4.8) 

where d'~{t) is the evolution in S^i-i with "+" boundary conditions, which starts from "— " at 
time t = 5ciL{log{5L)y . Thanks to monotonicity we can restart from "— " the evolution at time 
t = 5ciL{log{5L)Y and we used the hypothesis A^ to freeze all spins outside S^i-i to the value 
"+" in the time interval [5ciL(log(5L))'^, L^]. Via a trivial time translation, the upper bound 
(4.5) (applied with L replaced by 5L — i) shows that the right-hand side of (4.8) is smaller than 
c/ L? (c being the same constant which appears in (4.7)) which, together with (4.7), completes 
the inductive proof. 

Proposition 2 

D 
Proof of Proposition 3. Clearly, it is sufficient to prove the claim for L large enough. In the 
following, 5 will denote a small positive universal constant (conditions on its smallness will be 
specified later). Given a point x G M^, let (r := ||x||,^,(/)) be its spherical coordinates, where 
G [0, tt] is the polar angle and (p £ [0, 27r) is the azimuthal angle, with the convention that the 
half-plane {2: G M^ : z^^' > 0, z^"^' = 0} corresponds to (f) = and that ^ = identifies the positive 
z^^' axis. Note that the portion of dB^ contained in the octant C+ := {(0,(/)) G (0,7r/2)^} is a 
monotone surface, where: 

Definition 2. A smooth surface F C R^ is said to be monotone if for every x £ T all three 
components of the normal vector Ux at x are non-zero and have the same sign. 

It is convenient to introduce a few geometric definitions: 

Definition 3. For a given x G Sl, let 

(a) L be the infinite half-line which starts from the origin of M^ and goes through x 

(b) x' be the intersection of L and dBi 

(c) n be the plane perpendicular to L which meets L at distance L — (logL)^'^ from the origin 
(the exponent 3/2 is somewhat arbitrary, and it could be replaced by 1 + 5 for any 6 > Q) 

(d) T be the half-space not containing the origin and delimited by 11 

(e) M. be the spherical cap Bl n T. 

Also, decompose the boundary of Ai as dAi = Fi U F2, where Fi = dM n dB^ is the "curved 
portion of the boundary" and F2 = dAi Pi 11 is a disk. 
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Thanks to the discrete symmetries of the sphere Sl, it is clearly enough to prove (4.4) for 
points X £ Sl\ Sl~i such that 9 G [0, it/2] and (p G [0, vr/4]. Actually, we claim that it is enough 
to restrict to 9 £ [0, 7r/2 — 6] and (j) G [0, 7r/4], if 5 > is small enough. Indeed, if we interchange 
the role of the first and third coordinates of M^, the points with cp £ [0,-7r/4] and 9 close to -7r/2 
are mapped into points with 9 < 7r/4 and (j) small. 

It is convenient to distinguish three cases (see Figure 1): 



CaseB 



Case A 




Figure 1. The monotone octant C"*" of the surface of the sphere of radius L. 
By symmetry, we need to prove (4.4) only in the case where <j) < 7r/4 and x' 
(cf. Definition 3(b)) does not fall in the strip of width 6L, adjacent to the {x,y) 
plane. Case A corresponds to x' in Cl (the subset of the surface of the sphere 
delimited by the thick line). Case B corresponds to x' at geodesic distance at 
most \'LlogL from the north pole. Finally, case C corresponds to x' in the strip 
of width (1 + 6)\/2L(logL)^'^ to the left of Cl- For reasons of graphical clarity, 
proportions are not respected in the drawing. 



Case A: X £ Sl\ Sl-i is such that (6 

Cl ■■= {[0,n/2 - 5] X [0,tt/A\) n I {9,, 



£ Cl where 
logL ^ 



> 



L 



sin(^)>ii±^(logL)3/^ 



L 



(4.9) 



Remark that the lower bound on (j)sm{9) just means that, moving on OBl along a line of 
constant 9, the distance between x' (cf. Definition 3(b)) and the set of points in OBl where 
(/) = is at least (1 + <5)\/2L(logL)3/4. 

Lemma 2. Under the condition {9,(j)) £ Cl, one has: 

(a) The disk Fa has radius ^/2L{\ogLf/'^{l + o(l)) 

(b) Fi is contained in the monotone octant C~^ = (0,7r/2)^. 

As a consequence, both Fi and F2 are monotone surfaces, cf. Definition 2. 
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Proof of Lemma 2. Statement (a) and monotonicity of r2 require just elementary geo- 
metric considerations. As for statement (b), it is easy to see that the condition {6,(p) £ Cl 
imphes that the geodesic distance of x' along dB^ from the boundary of C"*" is larger than 
(1 + (5)v2^(logL)^'^(l + o(l)). The fact that Fi is contained (for L large enough) in C^ then 
just follows (thanks to the fact that S is strictly positive) from statement (a). 

Lemma 2 

□ 

We can now continue the proof of Proposition 3 under the condition {9,(/)) G Cl. Let 

be the collection of lattice sites contained in Ai (cf. Definition 3(e)). It is clear that the site 
X G Sl\ Sl-i under consideration belongs to Ai and that, thanks to Lemma 2(a), 



D := max 

y,z€M 



[d ((zW, z(2)), (y (1), y(2))) } < 2V2L{logLf/\l + o(l)). (4.10) 



By monotonicity, it is enough to prove (4.4) for a modified dynamics where only the spins 
y £ Ai evolve, starting from the "— " configuration, with boundary conditions given by r^j^ = + 
for y ^ Sl and i]y = ~ for y £ Sl\M. For lightness of notation, we still call such dynamics 
a-{t). Let 

y+:=5LnN^ V-:=V+\M, (4.11) 

and observe that V are positive monotone sets in the sense of Definition 1, as guaranteed by 
the monotonicity of the surfaces r^,F^ (cf. Lemma 2). Also, call Mt := {z G N^ : cr^(t) = — }: 
thanks to the discussion in Section 4.1, we have that Mt is a positive monotone set at all times 
t > 0, and trivially V~ C Mt C y+. To be coherent with the notations of Section 4.1, we have 
in the present case Vq = V^, A = V^ \ V~ = Ai; also, the law of Mt is just u^ °, with invariant 
measure Py±- Finally note that, from the definition of the spherical cap Ai, one has 

i7:=max{|z(3)-2/(3)| , z,y £ Ai and y^"^ = z^^^a = 1,2} < ^^°^^j^ < 0(6) (log if ^\4:. 12) 

(recall that we are working under the assumption that 9 < tt/2 — 5). Putting together Theorem 

4 with the estimates (4.10), (4.12), we get that the mixing time of the dynamics in A^ is 
0(L(logL)i3/2). 

We have the following key equilibrium estimate: 

Proposition 4. Recall that py± is the uniform distribution over all positive monotone sets 
V G S+ such that V~ C y C T/+. There exists c > such that for every L > 2 and whenever 

3z£M: d{z,V-) > j{logLf/^ and z £ V < - exp(-c(logL)3/2). (4.13) 



Pv± 



c 



In words, this result is saying that spins which are at distance of order (logL)^'^ away from 
the bottom of the spherical cap A^, where the "— " boundary conditions act, are "+" with 
overwhelming probability. It is crucial here that the estimate (4.13) (i.e. the value of c) does 
not depend on the spherical coordinates {6, (p) of the point x under consideration, i.e. on the 
slope of the plane 11 of Definition 3(c). 

Roughly speaking, (4.13) follows from recent works on the height fluctuations of random 
monotone interfaces associated to dimer coverings of the hexagonal lattice [9, 8, 10], but it 
requires some work to really prove the precise statement we need. The proof of Proposition 4 
is given in detail in Section 5. 
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Once we have Proposition 4, the proof of (4.4) proceeds via a standard argument which we 
simply sketch: Eq. (2.8), together with the fact that the mixing time of the dynamics in Ai is 
0(L(log L)^^' ^), impHes that the variation distance between py± and the law of a~ (t) is smaller 
than exp(— c'(logL)^'^) for all t > cL(logL)^. Since the dynamics undergoes 0{L^) updates 
during the time interval [cL(logL)^,L^], via a union bound and the equilibrium estimate (4.13) 
one gets that the probability in (4.4) is lower bounded by 

1 - L^" exp(-c(logL)3/2) > 1 _ cL"^ 

for L large, which is the desired bound. 

Remark 4. It is important to notice that exactly the same proof gives, for some c > and for 
all L > 2, 

inf P {a~{t) = + for all times t G [cL(logL)'=, L^]) > 1 - ^ (4.14) 

^e5L\'S^_(3/4)(l„gi)3/2: L 

(just look at the equilibrium estimate in Proposition 4). 

Case B: X £ Sl\ S'l-i is such that 9 < L'^/^ logL. 

Here, a rather rough argument suffices. Remark first of all that the vertical coordinate x'^-* 
of X belongs to the interval [L - (l/2)(logL)2(l + o(l)),L]. Call S := {y e Sl : y^^'^ > x^^^}. 
By monotonicity, to show (4.4) it is sufficient to prove that 

P(r+ > cL(logL)") < c/L^ (4.15) 

for a modified dynamics where only spins in S evolve, starting from the "— " configuration, 
with boundary conditions given by 77^ = + if z ^ 5^ and rjz = — ii z £ Sl\ S. Note that S 
is a discrete spherical cap, and that the boundary conditions just described are "— " below its 
base, and "+" elsewhere. Of course, r-f in (4.15) is understood to be the first time when all 
spins in S are "+" . 

We note first of all that S is contained in a parallelepiped Q whose base is a square of side 
£ := 2\/LlogL(l+o(l)) and whose height is /i := (l/2)(logL)^(l+o(l)). Next, by monotonicity 
we see that r+ is stochastically increased if we replace S with Q, with "— " boundary conditions 
below its base square, and "+" everywhere else. One can decompose Q into h horizontal squares 
of side i, stacked one on top of the other. If h were 1, we would simply have the evolution 
for (3 = +CX) of the two-dimensional Ising model with "+" boundary conditions in a square of 
side i (the "+" boundary conditions on the top face of Q would compensate exactly the "— " 
boundary conditions on the bottom face), and we know [7, Theorem 1.3 (a), case d = 2] that 
in this case one has P(r+ > ci"^) < exp(— 7^) for some positive 7, if c is large enough. Via 
a standard monotonicity argument (cf. [7, Proof of Theorem 1.3 (a), case d > 2] for details) 
this implies that, for the evolution in the parallelepiped Q, it is very unlikely that r_|_ exceeds 
h X cf: 

P(r+ > cL(logL)^) = P [r+ > cfh{logLf-'^{l/A + o(l))] < he^-^/l) (4.16) 

if c > 4, which is actually stronger than the estimate (4.15) we wished to prove. 



Case C: X £ Sl\ Sl-i is such that 






and 6 sin 9 G 



0.(i±^(,„gL)3/- 



(4.17) 



14 P. CAPUTO, F. MARTINELLI, F. SIMENHAUS, AND F.L. TONINELLI 



Definition 4. Let x he the point of dBi with angular coordinates 

{«,* + «:= fe,*+^2^) (418) 

where 7 > will be chosen later, and let x be a point ofTs^ of minimal distance from x. 

For every y G M^, let (f(y), 6{y), 4>{y)) be the spherical coordinates of the vector y + 3x, i.e., 
the spherical coordinates of y with respect to the point —2>x G Z^. 

Lemma 3. We have, for the point x £ Sl\ Sl^i under consideration, 

f{x) = 4L - ^(log(4L))3/2(l + 0(1)), (4.19) 

e{i + o{i)) < e{x) < e (4.20) 

^{x)sme{x) = (/.sine + ^^(log(4L))3/^(l + 0(1)). (4.21) 

We omit the proof of Lemma 3, since it requires only tiresome but elementary computations: 
one first writes down the Cartesian coordinates of x and x, then one works out the spherical 
coordinates of x + 3x; finally, the three statements are obtained by expanding these spherical 
coordinates for L large, using the fact that 0' <C 1. 

We can now conclude the proof of (4.4), case C. First of all, remark that, thanks to (4.19), 
one has for L large enough and letting Sr,a = ^'^ H -Br,a, with B^^a the ball of radius r centered 
at a, 

X G 5'4L _3s \ 5'4^_(3/4)(log(4L)3/2) -3S 

provided that 

¥ < \ (-^) 

Secondly, (4.21) implies that for L large enough 



(x) sin6'(x) > 

yjAL 



a±P(.ori4L))'/' 
provided that 

^>(1 + (5)^/2. (4.23) 

Note that conditions (4.22) and (4.23) are compatible if 5 is small enough, and choose a value 
for 7 which satisfies both. Finally, (4.20) shows that 6 < it/2 — S. 

All in all, thanks to Eq. (4.14) in Remark 4 (applied with L replaced by 4L, and modulo a 
trivial translation of the center of the sphere S/j^l from to — 3x), we have 

P (a^(t) = + for ah times t G [cL(logL)^L^]) > 1 - c/(4L)^ (4.24) 

with P the law of the dynamics in the sphere S4L-3X with "+" boundary conditions on 
dS4L,-3x, started from "— ". Since S4L-3X ^ Sl, by monotonicity (4.24) implies the same 
inequality when P is the law of the dynamics in Sl, always started from "— " and with bound- 
ary "+" conditions on OSl. The desired estimate (4.4) is proven. 

Proposition 3 

D 
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5. Estimates on height fluctuations of monotone interfaces 

Definition 5. A subset V QT? is said to be monotone if x ^ V implies y €z V whenever y G Z^ 
and y^"'> < x'"', a = 1, 2, 3. The collection of all monotone sets is denoted by S. 

Remark that this definition differs from that of positive monotone set (cf. Definition 1) only 
in that it is not required that y C N'^. Given a positive monotone set V G S'^, in the following 
we will always implicitly identify it with the monotone set V U (Z^ \ N^) G S. 

Definition 6. Given V £ T,, we associate to it a vertical height function, which we denote v. 
The function {vx}x&? ^* defined as 

v^ :=max{x(3) : (x(^\x(2),x(^)) G y} i/ x = (x^^^x^^^). (5.1) 

Observe that Vx takes values in Z U {— oo, +oo}, and that 

Vx < t-j/ify^"^ < a;("\a= 1,2. (5.2) 

We will denote V the set of all possible functions f : Z^ i— ;■ Z U {— oo, +00} which satisfy (5.2). 
As discussed in Section 5.2, one can identify V with 2?-^ x Z, where P-^ is the set of dimer 
coverings of the infinite honeycomb lattice Ti. 

5.1. Proof of Proposition 4. Recall the definition (4.11) of the monotone sets V^ in Section 4, 
and note that the corresponding height functions w^ G V coincide outside some domain {/ C N^ 
whose diameter is 0(vL(logiv)^'^) (cf. Lemma 2(a)). Note that U is just the projection on the 
plane (x,y) of the discrete spherical cap Ai. 

The estimate (4.13) is proven if we have, for some universal constant c > 0, 



Pv'k 



. {A\r^ n r^) < - exp(-c (log Lf^^) (5.3) 



for L > 2 where 



A := {ax. t/: „.>.- + 2^} (5.4) 

r := {f ~ < Vx for every x € U} (5-5) 

r := {v^ > Vx for every x G t/} (5-6) 

and py\ ^(•) (for some v £ V) denotes the uniform measure over the elements w of V such that 
Vx = Vx for X ^ U. (Since monotone sets and height functions are in one-to-one correspondence, 
we use the same notation p to denote the equilibrium uniform measure in both cases). Note 
that the event A is increasing with respect to the natural partial order in V, where we say that 
v < w (with v,w £ V) if Vx < Wx for every x G Z^. To lighten notations, given yl C Z^ we will 
write VA < wa if Vx < Wx for every x G ^. 

We need the following monotonicity property, which is an immediate consequence of Propo- 
sition 1: 

Lemma 4. One has 

Pwl^A-lO'U <wu < bu) < Py,i\^S-Wu <wu < b'u) (5.7) 

if a,a',b,b' ,w,w' G V are such that au < a'jj, bjj < b'jj and wjj^ < w'^jc. Moreover, p^| ^(•) 
depends only on the value of w on dU . 
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Applying Lemma 4, we get 

p,-i^^(^|rinr2)<p,-|^^(A|ri). (5.8) 

The key point is the fohowing result, whose proof is given in Section 5.2: 

Theorem 5. There exists a probability measure v on the elements tf G V which satisfies the 
following properties: 

z/(3x £dU ■.w^<v-)<eL (5.9) 

^{PM^SA))<eL (5.10) 

KPHac((r')'^))<eL, (5.11) 

where 

EL ■■= (1/c) exp(-c (log L)3/2) (5.12) 

and c is a universal positive constant. 

We are now in a position to prove Proposition 4. We have 

'^(pwl^A^l^^)) > J^ Pw\^AA^^)\VdU^^dU '^(Vdu^^du) (5.13) 

where we used (5.9) and Lemma 4. Therefore, 

p.-l^AA\r') < 2z.(p^|^^(A|ri)) (5.14) 

= 2u{p^\^^ {A\T'y, p^i^^ (pi) < 1/2) + 2u{p^\^^ {A\T'y, p^,^^ (F^) > 1/2) 

< 2u{p^\^^ {{TY) > 1/2) + ASL < MpMuc ((r')"^)) + 4^^ < 8£L 
where we used assumptions (5.10), (5.11), Markov's inequality and the obvious bound py,\ ^ [A\T^) < 

Eq. (4.13) 

p^l ^{A)/pyj\ ^(r ). Recalling the definition of El, this implies (4.13). D 

The idea behind Theorem 5 is that, while the height fluctuations of v under p for the Exed 
boundary conditions imposed by f~|c/c are hard to control, they are instead easily described 
if the boundary conditions are sampled from the infinite measure u described in Section 5.2. 
Such random boundary conditions are with high probability "higher" than the deterministic 
ones (see (5.11)) and an application of monotonicity (cf. the steps in (5.13)-(5.14)) concludes 
the argument. 

5.2. Proof of Theorem 5: Dimers coverings and height functions. We need first to recall 
some notions and results about the connection between monotone sets and dimer coverings of 
the infinite two-dimensional hexagonal lattice (cf. for instance [9, 8, 10, 11]). 

For X G M^, let '7riii(x) denote its orthogonal projection on the (111) plane {z G M^ : z'^'^' + 

T = U^ez37riii(z) (5.15) 

and note that T is a two-dimensional triangular lattice, which we consider as a graph by putting 
an edge between any two nearest neighbors (whose mutual distance is y^2/3). We also let the 
vector ei (resp. 62 and 63) be the vrm projection of the vector which joins (0,0,0) to (1,0,0) 
(resp. to (0, 1, 0) and to (0, 0, 1)): of course, the Cj have norm a/2/3. See Figure 2. 
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Let Ti be the dual lattice of T, again with edges between nearest neighbors. 7^ is a hexagonal 
lattice which we decompose as % = T-Lb U Ti-w-, where both T-L\y and T-Lb are translates of the 
triangular lattice T, such that all nearest neighbors of every vertex in T-Lw (resp. in T-Lb) belong 
to T-Lb (resp. to T-Lw)- Vertices in T-Lw (resp. in T-Lb) are said to be white (resp. black). 

Conventions^ We embed T and Ti in M^, with the convention that '7riii(0) G T is mapped 
to (0,0), that 6=3 is mapped to the vertical vector (0, y^2/3) and that ei is obtained by 63 by 
a counter-clockwise rotation of (2/3)7r. Also, given the two endpoints of the edge of Ti which 
crosses 63, we decide that the one which has negative horizontal coordinate (call it wq^o) belongs 
to 7iw^ while the other (call it bo,o) belongs to TIb- Observe that the vector bo,o — wo,o is 
proportional to 62 — 61. We will label a site v G Tiw (resp. in Tis) as Wx^y if u = wq^o + a^ei + ye2 
(resp. as hx,y if f = bo,o + xii + ^62). Edges of 7i which are perpendicular to ei (resp. to 62 or 
63) will be called edges of type "a" (resp. of type "6" or "c"). See Figure 2. 

Given a monotone set V £ Ti, we define the height function h := h{V) := {hx}x£T as follows: 

hx := max{z(^) : 7riii(z) = x} € Z. (5.16) 

zev 

We denote by W the set of all functions h : T >-^ Z such that there exists V G 'S with h = h(y). 
It is easy to see that there is a one-to-one mapping between elements of S and elements of V 
(cf. Definition 6), and between elements of V and elements of W. In other words, the functions 
{vx}x&i^ (cf. Definition 6) and {hx}xeT are two equivalent ways to describe the height of the 
set V with respect to the horizontal plane. 

A dimer covering M of the hexagonal lattice ^ is a subset of the edges of 7i covering each 
vertex of Ti exactly once. Note that each edge in M covers one black and one white vertex. 
To each height function (or monotone set) /i G W is uniquely associated a dimer covering, and 
conversely to a dimer covering one can associate uniquely a height function, provided that one 
fixes the height function at some arbitrary point x £ T (in other words, dimer coverings identify 
only gradients of the height function) . 

The construction of the dimer covering M given h{V) goes as follows (see also Figure 3). Let 
e be an edge of Ti, and let {x,y) be the edge of T which intersects e, with the convention that 
the vector x — y is +ei for some i G {1, 2, 3}. li x — y = ei or x — y = 62, then we put a dimer 
on the edge e if hx = hy — 1, and we do not put it if hx = hy (it is easy to see that these are the 
only two possibilities since h is a monotone interface). If x — t/ = 63, then we put a dimer on e 
if hx = hy and we do not put it if hx = hy + 1. Note that the asymmetry between the indices 
1, 2, 3 is due to the fact that we are computing heights with respect to the horizontal plane. A 
more symmetric choice (but less convenient for our purposes) would be to measure heights with 
respect to the (111) plane. 

Conversely, the construction of h given a dimer covering M goes as follows. First we define a 
Bux Lo, i.e. a function on oriented edges e of Ti, such that uj{e) = —uj{—e). If e is oriented from 
the white to the black vertex, then: 

• if e is of type "a" or of type "6", then a;(e) = if e ^ M and uj{e) = 1 otherwise 

• if e is of type "c", then uj{e) = if e G M and uj{e) = —1 otherwise. 

Next, we fix some arbitrary value /i^ G Z at some point x £ T- Finally, the difference /i„ — hx 
for t; G T is the total fiux of uj which crosses, from right to left, a path i of edges of T, which 



The conventions we adopt on the orientation of the axes and on the labehng of the edges via the symbols 
a, b, c do not coincide with those of [9, 8], but this is simply an irrelevant matter of convention. 
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Figure 2. A portion of the triangular lattice T and of the hexagonal lattice %. 



goes from x to i;. The fact that h^ — hx does not depend on the choice of the path i is due to 
the fact that the flux u has zero divergence, see [9, Sec. 2.2]. 

Given p = {pa,Pb,Pc) with Pa,Pb,Pc > and Pa + Pb + Pc = 1, take a triangle of perimeter 
1 whose angles are 9a := Trpa, Ob '■= T^Pb, ^c := t^Pc and let ka (resp. kf,, kc) be the length of the 
side opposite to 9a (resp. 9b, 9c). To every choice of p as above, one can associate a translation- 
invariant Gibbs measure ^p on dimer coverings of %. Translation-invariance means that, if A 
is a set of edges of T-L, then /Xp(A C M) = fip{T{A) C M), where T is a translation which maps 
Ti into itself. 

The precise statement is the following, whose different pieces were proved in [11, 21, 10]: 

Theorem 6. There exists a unique translation-invariant law ^p on dimer coverings, such that 
the probability that a given edge of type "a" (resp. of type "5", "c"J belongs to M is pa (resp. 
Pb,Pc) and such that, conditionally on the configuration Mx-= of the covering M outside a given 
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Figure 3. A graphically convenient way to visualize a monotone set V is to 
associate to every x G V a unit cube, centered at x — (1/2, 1/2, 1/2). In this way, 
a monotone set, when seen from the 111 direction, appears as a tiling of a portion 
of the plane with three types of lozenges (top drawing on the left). If we mark a 
segment along the longer diagonal of each lozenge (top drawing on the right), we 
obtain a dimer covering of a portion of the hexagonal lattice (bottom, left). To 
each vertex of a lozenge, i.e. to every vertex in a triangular lattice, one associates 
its vertical height (bottom, right). One can check in the drawing that the dimer 
covering and the height function are linked by the construction explained in the 
text. 



domain X CT-L, Hp is the uniform measure over all coverings Mx of X compatible with Mx^, 
i.e., such that Mx U Mx^ is a covering ofH (we refer to this property as "DLR property"). 

Explicitly, /Up is described as follows. Define the m,atrix K := {i^(b, w)}{beWs,we'Hvi/} '^^ 
follows: 
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• ifhis not a nearest neighbor ofw, then K(h,w) = 0; 

• K{h,w) = ka (resp. kh,kc) if the edge (b, wj is of type "a" (resp. of type "6", "c"j. 

Define also the matrix K~^ := i^~^(w, b)|„g^j^ beWg} o,s 

K'^i'^x,y,K',y') = K'^{'^Ofl,bx'-x,y'-y) (5.17) 

and 

(27rz)"= Jy kc + kaZ + kbW z w 

where the integral is taken over the two-dimensional torus T := {(z, ty) G C^ : \z\ = \w\ = 1}. 
Then, given a set of edges X = {(wi,bi), . . . , (wfc,bfc)} in %, one has 

ti^{X C M) = \\{K{h,,^i)\ det (K-i(w„b,))^<^ .<,. (5.19) 

Note that X is a weighted version of the adjacency matrix of "H. It is immediate to check 
that one has the relation KK~^ = I, which justifies the notation K^^. The infinite matrix K 
however does not admit a unique inverse, as discussed for instance in [8]. 

Remark 5. If we fix deterministicahy the value hx for some x £ T, then the function T 3 x ^^ 
A*p(/ix) — hx is linear in x — x, thanks to translation invariance of ^p . As a consequence, the 
set of points 

{z G M? : ^in(z) G T, z^^^^ = Mp(/i.,n(.))} (5-20) 

is contained in some plane lip C M^ which of course contains the point z G M^ such that 
^111 (^) = X and z^^' = hx- It is rather easy to check that the normal vector of the plane lip is 
parallel to p. In particular, the plane lip is monotone in the sense of Definition 2. 

Height fluctuations. Set for lightness of notation A = (l/4)(logL)^'^. 

Proposition 5. Let x £ T and fix hx to some deterministic value. There exists c > such that 
for every L > 2 one has 

fip (3x G T such that d{x, x) < L and \hx - ^p(/ix)| > A/2) < - exp(-c(logL)^/^) (5.21) 

uniformly in p (where d(-, ■) denotes the Euclidean distance in T). 

Proof of Proposition 5. By translation invariance of ^p, we can assume without loss of 
generality that x = and that we fixed hx = 0. If d(x, 0) < L, we can write x = mei + ne2 for 
some m,n £ {—L, —L + 1, . . . , L}. Therefore, via a union bound, to show (5.21) it is enough to 
prove that, for every |n| < L 

/"p(|/ine, - /"p(/i„ej| > A/4) < - exp(-c(log L)^/^) (5.22) 

for i = 1,2. We consider for instance the case n > and i = 1, the other cases being essentially 
identical. From the construction of the height function in Section 5.2, we know that 

ho - /i„ei = |Mn{(bi,o,wi,i),...(b„,o,w„,i)}| =: 7V„, (5.23) 

i.e., it is just the number of dimers in the set {(wi^O)bi^i), . . . (w„^0)b„^i)}. Thanks to the 
determinantal representation (5.19), one has [22] that Mn has the same law as a sum of n 
independent Bernoulli random variables Xi,i < n, whose parameters qi := P{Xi = 1) are the 
eigenvalues of the matrix A := {ka K^^(wi^i, bj,o)}i<i,j<n. 
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In general, it is not easy to solve explicitly the double integral (5.18) which defines K~^. 
However, there are special values of x,y for which K~^{wQfi,hx,y) takes an easy form. In 
particular, one checks that^, for n G Z \ {0}, 

i^-^(wo,o,b„,-0 = (-ir^^^^ (5.24) 

and that 

i^"'(wo,o,bo,-i) = ^. (5.25) 

As a consequence, the entries of the matrix A are given by 

^~ > A^-3) \ ^ \ (5.26) 

Va « = J 

and then 

n 

Var(AA„) = J];(?,(l-(?i) = Tr(^)-Tr(^2) (5^27) 

i=l 

= nao(l-ao)-a?(2n-2)-a^(2n-4)-...-2a^_-^. (5.28) 

We will show in a moment that 

Var(7V;) < clogn (5.29) 

for some c < oo, uniformly in p, which allows to conclude the proof of (5.22): via the exponential 
Tchebyshev inequality, 

A^p(AA„-^p(AA„)> A/4) = P j J](X,-E(X,))>A/4] (5.30) 

\i<n J 

i<n 

Using the inequality exp(a;) < 1 + x + x^ which holds for — 1 < x < 1 and the fact that n < L, 
one deduces 

P{^fn-^ip{^fn)>A/A) < e-^/^ll{l + qi{l-qi)) (5.32) 

i<n 
< g-A/4+Var(Ar„) ^ g-(l/16)(logL)3/2+clogL_ ^^^^^-^ 

Analogously, one obtains the same upper bound for P{J\fn — fJ-pi-f^n) < ~^/4) and (5.22) is 
proven. 

It remains only to prove the estimate (5.29). Essentially the proof can be found in [9, Sec. 
6.3], where however uniformity with respect to p was not discussed, so we repeat quickly the 
necessary steps here. Observe first of all that 

laA < , ,^ , i G Z (5.34) 



The formula (5.24) differs from the analogous one in [9, Sec. 6.3] by the factor (—1)", probably due to a 
typo there. In any case, the global sign is inessential for our computation, since (5.27) below depends only on the 
absolute value of ii'~^. 
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uniformly in p. Therefore, 

Var(AA„) 



< 



n 



n 



ao(l -ao) -2^, 



«o(l - oo) 



oo 



n-1 



+ 2 ^ i af + 2n ^ ( 



j=i 



+ c log n 



(5.35) 



(5.36) 



for some c independent of p. Finally, one observes that the quantity in square brackets is 



identically equal to zero. To see this, let f{x) 
then observes that 



|x|(7r — |x|) so that ao(l — ao) = f{9a)/TT. One 



Ck 



1 

2^ 



f{x)e"'''dx 



!+(-!)'' 
P 



if /c = 

if /cGZ\{0} 



(5.37) 



-ikOa 



and then a straightforward computation shows that the Fourier identity f{9a) = '^kei "^^ ^ 
is equivalent to ao(l — ao) = 2 XlieN '^1- 

Proposition 5 

n 

Proof of Theorem 5. Recall the setting of Section 5.1. Let u = {u^^',u^'^') G dU, let z := 
(li(^) , u(2\ w^ + [A/(2cos(6'))J) G Z^, and x := 7riii(z) G T. Choose p to be parallel to the 
normal to the plane 11 of Definition 3(c) and consider the infinite volume distribution of height 
functions induced by fip, with the normalization hx = v~ + [A/(2 cos(^))J . Note that the plane 
Hp defined in Remark 5 is obtained from 11 by translating it for a distance [A/(2 cos(^))J +0(1) 
in the positive vertical direction, so that d{Il, Up) = A/2 + 0(1). 

For a given realization of the height function h, we let Y be the finite set of points 



Y{h) := {z G Z' such that d{TTiii{z), x) < L and z^ ' = /i7riii(z)} ^ ' 
Thanks to Proposition 5, one has that 

^p(3z G Y{h) such that d{z,Ilp) > A/2) < - exp(-c(logL)^/^) 



(5.3^ 



(5.39) 



with c independent of p. 

Recall that to a height function /i G W there corresponds a unique element v = v{h) G V and 
let u be the law on V induced by /Up. Note that V(u(i)^„(2)) = w^ + [A/(2 cos(^))J deterministically. 
We show now that v satisfies the conditions (5.9)-(5.11). The point is that, as observed at the 
beginning of Section 5.1, the diameter of [/ U dU is o{L), so that the ttih projection of any 
point {x^^'^x^"^' ,Vi,^(\) ^(2)\), with {x^^'^x^"^') ^ UU dU , has distance from x smaller than L. As 
a consequence, given some v = v{h) G V, 



B{v) := {x G 
and, from (5.39), 



such that 



(x(i), x(2)) GUUdU and x^^^ = v^^w,^^)} C Y{h) (5.40) 



1 



ly {3z G B{v) such that d{z,Up) > A/2) < - exp(-c(logL) 



3/2 X 



(5.41) 



which immediately implies condition (5.9), since the distance between Up and 11 is A, and the 
graph of the function U U dU 3 x >-^ v~ is within distance 0(1) from the plane 11 (recall that 
v~ is the vertical height function of the monotone set V~ defined in (4.11) so that, in U, v~ 
is just the lattice approximation of the plane H). Conditions (5.10), (5.11) are also immediate 
from (5.41), once one realizes that i'{py\i/c{0)) = 1^(0) if O is an event which depends only on 
{'Vx}xeu- This is because, thanks to the DLR property of /Up (cf. Theorem 6), the measure 
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i'{-\{vx}x£U'') is the uniform measure over the elements of V which coincide with v outside U; 
in other words, z^(-|{fx}xet/=) is nothing but p^|^j(-). 

Theorem 5 

D 



6. On the mixing time of a dynamics of monotone sets 

In this section we prove Theorem 4. Recall the notation from Section 4.1, in particular the 
definition of A = V~^ \ V~ , of D (the diameter of the horizontal projection of A) and of H 
(the maximal vertical distance between two points in A with the same horizontal projection). 
Without loss of generality, we can assume that A is contained in some parallelepiped F C N^ 
defined by 

r = {z G N^ : 1 < z(i) < ai, 1 < z^^^ < as , /ii < z'-^'> < /is} , (6.1) 

with ai < D and /i2 = maxlz*-^', z G A}, hi = nim{z^'^' , z G A}. Recalling Definition 6, one can 
identify uniquely any monotone set V G S"^, such that V" C T^ C V^"*", via a vertical height 
function vi^,j.,y)^ indexed by pairs (x, y) of non-negative integers. Since the configuration of V out- 
side A is fixed and coincides with y+, what matters is the collection v := {i'(x,j/)}i<x<ai,i<?;<a2- 
This defines a plane partition in the box T (cf. also Definition 1), i.e. a collection of heights such 
that V(^x+i,y) < '^{x,y)-: V(^x,y+i) ^ "^{x^y)- ^s discussed in Section 5, the bijective correspondence 
between v and V is given by 

V{x,y) = max{z : (x, y, z) ^V} . 

Each pair (x, y) can be identified with a unit square in the plane z^^' = 0, in such a way that 
the center of this square is given by the point of M^ with coordinates {x — 2^y ~ 2'^)' ^^^s, 
we interpret vtx,y) &s the height of the column at (x, y). The minimal and maximal sets V~ , V~^ 
correspond to minimal and maximal column heights, denoted v" and v~^ respectively (the same 
notation was used in Section 5.1). We define Q, = Q{v ) as the set of all plane partitions v in the 
box r such that v~ < v < v^ . Thus, fi is in one-to-one correspondence with the set {V £ S"*" 
such that V~ C F C V~^}, and the measure p = py± of Section 4.1 becomes now the uniform 
probability measure on Q (we still call it p). 

6.1. Column dynamics. The first step in the proof of Theorem 4 consists in establishing 
a mixing time upper bound (Lemma 5 below) for a Markov chain that involves equilibration 
of full columns at each move. The key idea here borrows from Wilson's analysis [25] of the 
Luby- Randall-Sinclair Markov chain for lozenge tilings [15]. The second step (see Section 6.2) 
is to show that Lemma 5 implies the upper bound on the mixing time of the "single spin-flip" 
dynamics v^ " which is under consideration in Theorem 4. A similar strategy was used in [17] in 
the simpler context of the (1 + l)-dimensional SOS model. 

A column (x, y) is said to be even/odd if (x — y) is even/odd. Consider the continuous time 
Markov chain with state space ri, where at each arrival time of a Poisson process with parameter 
1 we flip a fair binary coin; if the coin is (resp. 1) we update simultaneously all even (resp. odd) 
column heights vt^.y) with a sample from the distribution p conditioned on the current value 
of the height of the odd (resp. even) columns. Let Pt{v, •) denote the distribution at time t of 
such Markov process when the starting configuration is f G il. Note that the kernel Pt = Pt{-,-) 
satisfies 

Pt = e*^ , (6.2) 
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where the infinitesimal generator Q acts on functions / : fi i— t- M by 

[^/](^) = ^/>[/|odd](t;) + ^/>[/|even](t;)-/(^), (6.3) 

where we use the notation p[f \ odd](f) (resp. p[f \ even](ti)) for the expectation with respect 
to p conditioned on {t'(a;,y) , ix,y) odd} (resp. {'V(x,j/) > ix,y) even}). Note that samphng from 
p[- \ odd ] amounts to pick uniformly at random a new configuration of even column heights that 
is compatible with the current odd column heights (compatibility here means that the configura- 
tion of all column heights is then a plane partition v that satisfies v" < u < f "'"). It is important 
to remark that the probability measure p[- \ odd] is a product of single column probability mea- 
sures, i.e. conditionally on {v(^x,y) ; {x, y) odd}, all even columns become independent. The same 
remarks apply to p[- |even]. Clearly, p is the reversible invariant distribution of our Markov 
chain. 

Lemma 5. If H < D, there exists a constant c > such that 

sup \\Pt{v, •) - pII < cD^ exp ( - -^) . (6.4) 

Proof. We use the well known "lattice path representation" of a plane partition; see e.g. [15, 25]. 
Namely, any plane partition w G fi can be seen as a collection of /12 — /ii + 1 non-intersecting 
lattice paths 0'*^ = (j)^^\v), i = /ii, . . . , /12, each of length 2D [hi, /i2 are the same integers that 
appear in (6.1)) which satisfy 

4'^ = <P?D = J > ^li - 4^^ e {-i> +1} > 4'^ < ^i^"^ (6.5) 

for all J = /ii, . . . , /i2 , X = 0, . . . , 2D. The polymer j describes the j^^ level set {z G F n A : 
z = i}- For the precise construction of the paths, we refer to [25, Section 5], [15, Section 2.1] 
(see also Figure 4 for a graphical construction). 

In the plane partition-to-lattice path mapping, inequalities are reversed, i.e. if vi < V2 then 
(p{vi) > 0(^2) (cf. Figure 4). We let 0+ := ((^"'"'^*^)i and 4>~ := ((/>"' *^*'')j denote the lattice paths 
corresponding to the maximal and minimal plane partitions v~^,v~, so that (/>"*" < (f>~ . Then, 
the condition v~ < v < v'^ gives 

cPt'^^^ < 4J^ < rj^^ , (6.6) 

for all J = /ii, . . . , /i2 , X = 0, . . . , 2D or, more compactly, (/)~^ < </* ^ </*"• 

Let Q denote the set of configurations (j) = {4'x )x,j of integer heights (px G ^, J = ^1, • • • , /12 
and X = 0, . . . , 2D satisfying the constraints (6.5), (6.6). The construction above establishes a 
one-to-one correspondence between the set f2 of plane partitions v satisfying v~ < v < v~^ and 
the set ri. The image of the measure p is the uniform probability distribution on $7 (which we 
call again p with some abuse of language). Moreover, it is not hard to see that the image of the 
Markov process with "full column moves" under this map coincides with the continuous time 
Markov process with state space ft obtained as follows: at each arrival time of a Poisson process 
with parameter 1 we flip a fair binary coin; if the coin is (resp. 1) we update simultaneously 

all 01 , i = /ii,...,/i2, for X even (resp. odd) with a sample from the uniform distribution 
conditioned on the current values {4>x , i = hi, . . . , /12, x odd} (resp. even). With a slight abuse 
of notation we call again Pti-, ■) the kernel of the Markov process on lattice path configurations, 
and its invariant measure is of course the uniform measure p. Moreover, we write again (cf. 
(6.3)) 

g = ^(p[.|odd]-l) + ^(p[-|even]-l). 
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Figure 4. In the left drawing, the "column" or "plane partition" representation 
of a monotone set V G S"*". In the right drawing, the corresponding "lattice 
path" representation. The lattice path (j)^^' is essentially the level set at height j 
of the column representation of V, seen from the (111) direction. It is clear that, 
whenever we add a cube to the plane partition on the left, one lattice path on 
the right will be lowered. This shows that vi < V2 is reversed into 0(ui) > 4'{'^2)- 



for the generator of P(, where I denotes the identity operator. 
Next, we turn to Wilson's coupling argument. Define <I> : il i— )■ 



by 



2D~l 



$(</)) 



j=hi 



'^^ 



g{x) 



sm 



x=l 



/ vrxx 
\2D) 



To compute the action of Q on <1>, fix some even < x < 21?, and observe that if the constraints 



-,(i) 



< 



^(i) 



< 



-,(i) 



were absent, one would have 



h2 



odd 



3=hi 



hi 
j=hi 



e+ 



ei) 



(6.7) 



(J) 



Note that the sum is important in (6.7) since the identity has no reason to hold for a single j. 
Now, the constraint (6.6) is felt at x in the path j iff either A-{j,x) := {(/>; 
'^SJil or A+{j,x) :={(^: 



^.•^^■)< 



e^ 



U) 



l>p^'' > 



^(i) 



cannot move and 



(6.7) with —1 since 

have to compensate (6.7) with +1 since 
conclusion. 



J_i = (t>x^i}- In the first case we have to compensate 



dj) 



^X+lJ 



1. In the second case we 



cannot move and 



kiJ) 



e+ 






1. In 



h2 



j=hi 



I-') I odd 



. h2 
j=hi 



e 



h2 



1 + 4+i) + E (lA+a,.)('A) - IaA3,x){^)) ■ (6-8) 
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^h2 Aj) 



^h2 Aj) 



Clearly, for any x even one has p Xlii/n 'Px \ even = Yl,j=hi ^^ ■ Moreover, exactly the same 
equations hold when x is odd, provided we change the conditioning from odd to even. Therefore, 
using the notation {A(p){x) = ^{^{x — 1) + (p{x + 1)) — (p{x) for the discrete Laplacian, (6.8) 
together with its analogue for odd x imply that 

2D-1 h2 2D-1 /i2 

x=l j=hi x=l j=hi 

Summing by parts and using Ag = —ko g, where k^ = 1 — cos(7r/(2L')), one has 

2D-1 h2 

[Q^m = -'^H4>)+ E 5W Y.iU^{-j,x){<P)-U.M{4>))- (6.9) 

X=l j = hl 

Next, let ipAt) denote the state of the Markov process at time t with initial condition ^ G il 
at time 0. When ^ = (j)~^ (minimal state in terms of lattice paths) or ^ = (/)~ (maximal state), 
we simply write (p'^^t) or <p^{t). Define u{t) = E[$((/>~(t)) — ^{(/)~^{t))] > 0, where E denotes 
expectation with respect to the global monotone coupling of the lattice path Markov process 
with kernel Pt (cf. Section 2.1 and [25]). From (6.2) we infer 

|u(t) = -^n(t)+^(t), (6.10) 

with 

2D-1 h 

i^{t) = Y, 9{^)Y.^ [{u^ij,x0^{t)) - iA_o-,x)(0-(i))) + {wu,x){r{t)) - iA+a,.)(0+(t)))] 

x=l j=l 

By monotonicity of the coupling it is immediate to see that ip^t) < 0. Therefore (6.10) implies 
u{t) < u{0)e^^~ . Note that u(0) can be upper bounded by the volume enclosed between 
the minimal and maximal plane partition, i.e. u(0) < \V^ \ V^\ < D^H < Z)^ (recall that in 
Theorem 4 we are assuming H < D). It follows that u{t) < D^e~^~*. 

To finish the proof it suffices to observe that by monotonicity ^{(/)^{t)) — ^{(p'^{t)) > and 
0+(t) / <p-(t) iff $(</>"(*)) - $(</>+ (t)) > 2sin(7r/2D), so that by Markov's inequality 

p(</,+ (t)/</,-(t)) < ^^^L< ^^iL^e-T*. 

VS" V ; 7- s" \ jj - 2sm{TT/2D) " 2sin(7r/2Z)) 

We can now bound the total variation distance by the probability of no coupling up to time t, 
and using monotonicity and the bound above this gives, for any initial state v 

\\Pt{v,-)-p\\ < TT^^^TTT^e-^*, 
2sm(7r/2L') 

which is easily seen to imply the desired estimate. 

Lemma 5 

□ 

6.2. Proof of Theorem 4. The proof of Theorem 4 uses the bound of Lemma 5 together with a 
comparison argument that allows us to translate the mixing time upper bound of the "column dy- 
namics" to an upper bound for the mixing time of the dynamics Vf- ° defined in Section 4.1. In the 
plane partition language, the dynamics V^ ° will be called v^'-', with law f^° and equilibrium dis- 
tribution as usual denoted by p: it has initial condition u o G fi (vq being the plane partition cor- 
responding to the monotone set Vq)) to each column (x,y) is associated an independent Poisson 
clock with parameter 1 and the evolution proceeds by local updates. By this we mean that when 
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the clock of the column (x, y) rings, one replaces 'V(x,j/) by iaax{v(^x,y) ~ 1) "^(x+i.j/)) 'y(x,j/+i)) ^f^: v)^ 
(with probability 1/2) or by min{t;(^^j/) + l,V(^^_i^y),V(^^^y_i),v'^^^y^} (with probability 1/2). In 
other words, the chosen column performs a simple symmetric random walk step, except that 
jumps which violate the plane partition constraints Vf^^^y-^ > max{t;(j._|_x^y), ^(i^^y+i)}, 'U(x,j/) ^ 
min{f(2,_i ,^),f(2, j^^x)}; oi" the overall constraint v7 -, < ^(3,,^) < vt ., are rejected. We need to 
prove that the mixing time of this chain is O {D^H^ (log D)^). 

We start with a simple observation that allows us to reduce to the case of maximal {vq = v~^) 
and minimal {vq = v~) initial conditions. 

Lemma 6. For any t > and any wqi^o ^ ^■ 

II, /""o _ ,/o|| < r)3 II D+ _ i;-|| 
v' 

Proof. Let P denote a monotone coupling of t'^" , t'i " • Then, using the fact that each column 
(x, v) has a minimal height v7 . and a maximal height V/ . such that vt •, — v7 . < H < D 

V 'W/ to {x,y) to {x,y) {x,y) {x,y) — — 

we have 

<E E [nK^)(.,.)>/^)-nK")(.,.)>/^)]<^'ii^r-'^rii- 



Lemma 6 
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Thanks to Lemma 6, to prove Theorem 4 it is sufficient to show that 



for some t = O {D'^H'^ (log D)'^). Let us prove the statement for z/^ , the argument for u^ being 
identical. Consider i.i.d. Bernoulli(2) random variables ( = {(^1,(2, ■■ ■) and let T^ = iT, with 
T = ci H^ logD and ^ G N U {0}, denote a partition of the time axis. Furthermore, consider 
i.i.d. Poisson processes with parameter 1 at each column {x,y), and let {Si{x,y)}i(^^y^ denote 

the corresponding collection of arrival times. Call {Si{x,y)}i(,j.^y\ the collection of arrival times 
obtained from {Si{x,y)}i^(^x,y) by deleting (or "censoring") all arrivals Si{x,y) such that, for 
some j £ N, Tj_i < Si{x,y) < Tj and {x,y) has the opposite parity as (j (e.g. {x,y) is odd and 

Cj = 0). 

By construction, if we start from the configuration v^ at time and perform local updates 
using all the marks {Si{x, y)}i^(x,y) up to time t we obtain the distribution i^^ . Let us call uF^ 
the distribution of the "censored" dynamics obtained in the same way but only using the marks 
{Si{x,y)}i(x,y): for ^ fixed Bernoulli sequence (. From the "censoring inequality" of Peres and 
Winkler [20], [18, Th. 2.5] it follows that f^ is stochastically dominated by z^A/, for any t' < t. 

Moreover, setting V^" := lE^^'/'^ , where E^ denotes expectation over the random sequence ("> by 
linearity of the expectation one sees that 



't 



^ l^^f for any t' < t. (6.12) 

Let us fix to = C2sT, where s = c^D^ log D and 02,03 are constants to be taken sufficiently 
large. Let {N{s),s > 0} denote a Poisson process with parameter 1, and write P and E for 
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the associated probability and expectation. Monotonicity implies that the event A = {v £ Q : 
uV {v) > p{v)} is increasing [18, Th. 2.5]. Therefore, the censoring inequality (6.12) gives 

Il<-Pll = vl^{A)-p{A)<n^^^\^)T{A)\N{s)T<t,)-p{A) (6.13) 



P(iV(s)T<to) ^¥{N{s)T<tQ) 



^ \\^^n\s)t-p\\ ¥{N{s)T>tQ) 

- l-¥{N{s)T>tQ) ¥{N{s)T<tQ)' 

Taking C2 large in the definition of to above and using standard estimates for the Poisson random 
variable we can make ¥{N{s)T > to) smaller than D~'^. Therefore, we see that thanks to Lemma 
5 and (6.13), if C3 in the definition of s is chosen large the claim (6.11) is a consequence of 

m^}^UT-Psiv^r)\\<^,, (6.14) 

where Ps{-,-) is the kernel defined in Section 6.1, which involves full column equilibrations 
(cf. (6.2)). To prove (6.14) we need to compare column equilibration moves with local up- 
dates. Let us use the notation pq{v,-) and pi{v,-) for the probability kernels associated to 
/9[-|odd](f) and p[-\even]{v) respectively, see (6.3). That is, for any / : il 1— t- M one has 
for instance /9[/|odd](f) = J2v'enPo{'>JTv')f{v'). Define P^^ = [p^-^ ' " P(;„]{v,-), and note that 
Q'^ := 2~" X^rejo i|" ^Cn ^^ nothing else but a discrete time version of the kernel Psiv, •), i.e. 
Ps{v, •) = KQ^^,-. where E and A^(-) are as above. Since s = C3 D^ log-D, excluding an event of 
probability 0{D^'^) for some large p > we can assume that N(s) = n for some n < 2s. Recall 
that V}'^ denotes the law u^^ conditioned on the binary sequence (. The previous remarks 
imply that it will be sufficient to prove the upper bound 

sup W^cU - PtnW < TTm ^ «<2s. (6.15) 

Observe that, for the censored dynamics i^^ ^ , in the time interval Tj^i < t <Tj, all columns with 
the same parity as Cj are independently updated by local moves, i.e. on columns of that parity we 
have independent continuous-time simple symmetric random walks in segments (determined by 
the columns of opposite parity) of length bounded by H. It is standard that for some constant 
c > the mixing time on each column is bounded by c H^ and therefore after a time T we have 
the bound, for any d G {0, 1}, uniformly in the starting configuration v £ Q: 



Kri 



PcAv,-)\\<cD^ew(^--^y (6.16) 



For general n G N, by recursive coupling of the distributions involved, using (6.16) at each step, 
one has 

IFc;nT(•)-kl•••/>cJ(^'">•)ll<l-(l-cI)2exp(--^))^ (6.17) 

Since T = ci H^ logD and n < 2s = 0{D^ logD), we see that the desired estimate (6.15) follows 
for a suitable choice of ci. 

Theorem 4 

D 



29 



7. The mixing time in dimension d = 3 (lower bound) 
Here we prove the bound 

P(r+<LV(c log L))<c/L (7.1) 

for a suitable constant c. 

Let Cl be the cube {1, . . . , L/2}^ (we assume for definiteness that L is even) with boundary 
conditions rjx = + ii x ^ OCl with min(x^^'',x'^'*, x*-^^) = and rjx = — otherwise. In other 
words, the boundary conditions rj are "+" at the three faces of OCl which meet at the origin 
of Z^, and "— " at the other three. Denote by {s~{t)}t the zero-temperature Glauber evolution 
started from the "— " configuration (in order not to confuse it with the evolution cr~{t) in A^) 
and by P^ its law. One has 

Proposition 6. There exists c > such that 

( L"^ \ c 

P!1 I there exists t < — such that sT, r /r,r\(t) = + I < — . (7.2) 

<-i \ clogL (1,-L/2,1)V ' J - L "- ' 

Proof of Eq. (7.1), assuming Proposition 6. Since the set {x G Cl '■ s^it) = +} is a monotone 
set (cf. Definition 1) at all times, the event that Sq^^iiIO = ~ implies the event Sy (t) = — 

for all y G Cl such that y^^' = L/2. Therefore, (7.2) implies (using also symmetry among the 
three coordinate axes) 

Pj^ [3 t < — -,x G Cl with max{x^^\x^'^\x^^'^) = L/2 such that s~{t) = +] < y.(7.3) 

The cube Al = {1, . . . , L}^ can be seen as the union of eight disjoint sub-cubes C*^*-* of side L/2, 
with C'^' = Cl/2 while C^^\i = 2, . . . , 8 are suitable translations of Cl- Let as usual P denote 
the law of the Glauber evolution inside A^, with "+" b.c, started from "— " and let P' be law 
of the evolution, again started from "— " , where the spin configuration inside each C^^' evolves 
independently for different i, with b.c. given hy rjx = + for x G dC^"^' n dh-L and i]x = — for 
X G 5C'*' n A^,. It is clear that, until the random time 

ti :=inf{t > : 3i G {l,...,8},x G SC^*^ n A^ such that cr~ (t) = +} , (7.4) 

the two evolutions can be perfectly coupled, and that ti < t+. On the other hand, thanks to 
(7.3) and to the symmetry among the various cubes C'*^ (up to suitable translations of the origin 
and reflections of the coordinate axes) one sees that 



^ ^ 


1 < 8 X |. (7.5) 


clogLy' 




Eq. (7.1), given Prop. 6 

D 



The main ingredient in the proof of Proposition 6 is a "column dynamics" s^ (t) (analogous to 
the one used in Section 6.1) defined as follows. Start from the "all minus" configuration in Cl, 
s~(0) = — . Assign to each x G {1, . . . , L/2}'^ an i.i.d. Poisson clock of rate 2; when the clock 
labeled x = {x^^',x^'^') rings, we assign a new value to the collection of spins at sites z G Cl 
with {z^^),z^'^') = {x^^^x^"^'), by sampling it from the equilibrium distribution conditioned on 
the present value of all the other spins. In other words, we set to equilibrium the column of 
horizontal coordinates (x^^\x^^^), conditionally on the value of the neighboring columns. For 
every t we have the stochastic domination 

s'{t) :<r{t). (7.6) 
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Indeed, for /c G N let s~' (t) be the following dynamics. Set s~' (0) = — and, when the Poisson 
clock of the column labeled x rings, repeat k times the following procedure: 

• flip a fair binary coin; 

• if the coin gives "head" , then make a heat bath update at each of the sites of the column 
X with vertical coordinate belonging to 2N, one by one, starting from the bottom site; 

• if instead the coin gives "tail" , then make a heat bath update at each of the sites of the 
column X with vertical coordinate belonging to 2N + 1, one by one, starting from the 
bottom site. 

Here, "making a heat-bath update" at a site z means updating az according to the equilibrium 
conditioned on the value of the spins outside z. It is immediate to realize that the process 
{s''^{t)}t has the same law as {s~{t)}t, and that the law of {s~'''{t)}t converges to that of 
{s^{t)}t for k —7- oo. Also, the stochastic domination s^' (t) ^ s^'^^{t) is an immediate 
consequence of the Peres- Winkler censoring inequality [20], [18, Th. 2.5]. 
Call P the law of the column dynamics {s^{t)}t. One has 

Proposition 7. Fix c > 0. For L sufficiently large and t < L^/(c log L) one has 

^i~^li,L/2,i)it) = +)<<^'L'e-^y^'''^ (7.7) 

for some finite constant c' . 

Proof of Proposition 6, given Proposition 7. Let 

so that, thanks to Proposition 7 and to the stochastic domination (7.6) one has 

Bl^{H) < c' L^-^'/'^^l (7.9) 

Note that the desired bound (7.2) can be rewritten as P^ {H > 0) < c/L. One has 

Pl,^{H>0) = Pl^{H>L-''^^^) + Pl^{d<H<L-'/^^^) (7.10) 

where in the first term we applied Markov's inequality. Choosing c sufficiently large, one can 
make both terms in the last expression smaller than c/L. Indeed, in order that < H < L~'^'^'^^, 
there must be two times s, t with < t — s < L^^' ^^® such that the Poisson clock associated 
to the site (l,L/2, 1) rings both at times s and t. Via a simple union bound, and using the 
exponential form of the law of the intervals between two successive rings, one sees that the 
probability of such event is 0{L'^~'^'^'^^). 

Prop. 6, given Prop. 7 

D 
Proof of Proposition 7. Since the set {x £ Cl '■ s~{t) = +} is a monotone subset of Cl 
at all times, one can identify it (recall Section 6.1 and Fig. 4) with the set of L/2 paths 

{(j) {t)}j=i,...,L/2 of length L + 1, where the j*'' path is the collection {(px {t)}x=o,...,L and 
the following relations are satisfied: 



it) = <p^\t) = j , 4%it) - #(t) G {-1, +1} , #(t) < 4^-+i)(t) . (7.11) 
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At time t = one has {x G Cl '■ s^ (0) = +} = and therefore (px {0) = j + x if x < L/2 and 
(j)i^\o) =j + L-xifx> L/2. Defining 

2 ^/' r 1 

hAt):=jY.[^i'Ht)-j\, (7.12) 



u{t, x) := E(/ia;(t/2)) and reasoning hke in Section 6.1 (see also [25, Sec. 5]), one sees that u{t, x) 
satisfies the discrete heat equation 

iuit,x) = iAu)it,x) := »ft.^+i)+"(^^^-i)-2^(*.^) ^ x = l,...,L-l .^_^^. 

u{0,x) =xl^<L/2 + {L-x)lx>L/2 

with Dirichlet boundary conditions u{t, 0) = u{t, L) = 0. On the other hand, from the represen- 
tation of monotone sets as collections of paths, one sees that 

P(VW)(*) = +) = ^^a^).m-i=-nl=^[ - ^"%^"^ " I (7.14) 




< 

Therefore, it is sufficient to prove the heat-equation estimate 

l-u(t,L-l)<c'Le-^'/(32t)^ (7_;L5) 

uniformly for t < 2L^/(clog-L) and L large. While (7.15) can be obtained directly via Fourier 
analysis, we give a simple and more probabilistic argument. First it is well known (cf. for 
instance [7]) that, for x = 1, . . . , L, the quantity [1 -|- n(t, x) — u{t, x — l)]/2 coincides with the 
probability that there is a particle at site x at time t, for a symmetric simple exclusion process on 
{1, . . . , -L} with initial condition at time zero such that sites x < L/2 are occupied by a particle, 
while sites {L/2 + 1, . . . , L} are empty (each particle attempts with rate one to jump to one of 
its two neighboring sites with equal probability 1/2 and the jump is rejected if either the site 
is already occupied or if it lies outside {1, . . . , L}). In particular (recall that u{t, L) = 0) one 
has that [1 — u{t,L — l)]/2 is the probability that there is a particle at L at time t. By duality 
(cf. [13, Section II.3]), this equals the probability that a continuous-time simple random walk 
of rate 1 on {1, . . . , L}, started from site L, is in {1, . . . , L/2} at time t. The bound (7.15) then 
follows from standard random walk estimates: if Px' is the law of the continuous-time simple 
random walk X^ on {n G N : a — 1 < n < 6 + 1} started from x, one has 

Pl''{Xt < §) < P^'^iXt <\)< P.^ti^i < §) < P,'-^'"-{^s < t : \Xs\ > ^) (7.16) 

and the latter expression is easily seen (e.g. using the local central limit theorem) to be upper 
bounded by the r.h.s. of (7.15). 

Prop. 7 

□ 

8. Proof of Theorem 2 
As discussed in Remark 2, we only have to prove that, for the /? = +oo dynamics, 

P (r+ < coL^) < exp(-7L) (8.1) 
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for suitable positive constants cq, 7. Also, thanks to monotonicity, it is enough to prove this fact 
for the dynamics o"~(t) in the domain 

A := Al n {x G Z^ : \x^^')\ + \x^^'>\ < L + 1}, (8.2) 

with "+" boundary conditions on dA. The advantage of looking at the dynamics in A instead 
of Al will be apparent in the proof of Theorem 7. 

We need a certain number of geometric definitions: 

Definition 7. Given a G $7^, let 
(a) 

M{a):={x: a^ = -1}; 

(h) 

T{a) := \Jx(zM{a)Bx-, 
where Bx is the unit square of side 1 centered at x, with sides parallel to the coordinate axes. 
The boundary ofT{a) is denoted as dr{a) and its geometric length as \dT{a)\; 

(c) for i = 1,2 

«^L(^) := max{a;» : x = (x^^^x^^)) g dT{a)} 
and similarly 

u2„(a) := min{x« : x = (x^^^x^^^) G dr{a)}; 

(d) p{cr) be the configuration in $7^^ obtained by flipping every "— " spin in A such that a strict 
majority of its nearest neighbors are "+ ", and repeating the same operation as long as such 
a site exists. We call p{-) the "majority transformation"; 

(e) D:={xel? : |x(i)| + |x(2)| < (9/10)L} C A; 

(f) Q be the subset of Vi-^ defined by 

Q := {cr : i9r(cj) is a simple curve (i.e. without loops), (8-3) 

2 
|ar(a)| = 2^(^« Ja) - n2„(a)),Z) C M{a) andp{a) = a}. (8.4) 



1=1 



2 c.XO {„\ ,,(«) 



Note that the constraint |9r((T)| = 2 Yli=ii''^"iax{o') — Ummi'^)) ™ ^he definition of G is simply 
the requirement that a minimizes the Hamiltonian H^{-), given the values {urnax{cr),Uj^^j^{a)}i=i^2- 

Next we introduce a modified dynamics {o"^(t)}t>o on A , with initial condition o'^(O) = ^, 
as follows. Let r|, := inf{t > : M{a^{t)) ^ {D U dD)}, and set d^{t) := CT«(r|,) for t > t%. 
Whenever the Poisson clock labeled x rings at a time t < tJ~^, first refresh the current value of 
the spin at x according to the distribution vr(-|o"y = ay{t),y 7^ x), and then apply the majority 
transformation p{-) of Definition 7(d) to the configuration thus obtained. 

We call C the generator of such dynamics and fll its law at time t. When the initial condition 
is ^ = — , it is immediate to see that, by the monotonicity of the usual Glauber dynamics, one 
has 

{M{a~{t)) ^ (DUdD)} =^ {M{d-{t)) ^ (DUdD)} (8.5) 

for every t > 0, so that (8.1) is proven if we show that 

P(r^ < coL^) < e-'^^ (8.6) 

for some suitable co, where P is the law of the process {5'~(i)}t>o- 
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The advantage of the modified dynamics a (t) with respect to the usual Glauber dynamics 
is that it belongs to the "good set" Q for all times: 

Theorem 7. For every t > one has that a^it) G Q. Moreover, there exists a deterministic 
positive constant tp such that 

|M(a-(0))|-|M(a-(r^))|>VL'. (8.7) 

Proof of Theorem 7. For this, we need some additional notions: 

Definition 8. Given a £ Q and x £ A, we say that x is a "flippable site" if two neighbors of 
X, at mutual distance y/2, are "+" and the other two neighbors are "—". If x is flippable, recall 
that a^^> is the configuration obtained by flipping ax to —ax- Then, for x flippable we say that 

• X is a "mountain" if x £ M{a) and a^^' G Q; 

• X is a "vertex" if x £ M{a) and p{a^^') / a^^'; 

• X is a "valley" if x £ dM{a). 

Remark 6. Observe that, if x is a "vertex", then exactly one of its neighbors, call it y, is also 
a vertex. It is immediate to see that M{p{a^^>)) = M{a) \ {x,y}, cf. Figure 5. 

We prove that a~ (t) £ Q foi t < t^ by induction: clearly the statement is true at time zero, 
and we show that if it is true until some time s then it remains true after the next update. 

The first observation is that, if x •= '^~('^) ^ ^i then nothing happens in the evolution 
until the Poisson clock of a flippable site x rings (the occurrence of sites with three neighbors 
of opposite sign, or sites with exactly two neighbors of opposite sign at mutual distance 2 is 
forbidden by the condition x ^ 0)- When such a ring happens, we have three possibilities: 

(1) x is a valley. Then, with probability 1/2 the configuration remains unchanged, and with 
probability 1/2 we change x to x S cf. Definition 8. We do not need to apply the 
transformation p{-), since one sees easily that p{x ) = X ■ Also, it is easy to see that 
|r(x)| = |r(x'^'')|, and that T[x ) remains a simple curve, see Figure 6. 

(2) X is a vertex. With probability 1/2 the configuration remains unchanged, and with 
probability 1/2 we change x to X i then, the application of the transformation p{-) has 
the effect of flipping also the vertex neighbor of x, cf. Remark 6. Altogether, F remains a 

simple curve; its length decreases by 2, but also does the sum 2 Yli=ii'^rnax{x)~'Umin(x))j 
see Figure 6. 

(3) X is a mountain. This case is more subtle, since it is not obvious apriori that T{x^^') 
is a simple curve. Just to flx ideas, assume that the "+" neighbors of x in x are 
x + (0, 1), x + (1, 0). If T{x ) were not a simple curve, it would mean that Xa;-(i,i) = +• 
Since x ^ G, one has that the set 

Y{x) := A n ({x - (1, 1) - {i,j),i,j > 0} U {x + {i,j),i,j >0,i + j> 0}) (8.8) 

belongs to A \ M(x) (so in particular it has no intersection with D), otherwise the 
condition |9F(x)| = 2 ^^^^(^^^^^^(x) - Ujninix)) would be violated. However, by the 
definition of the domains A and D, there exists no site x G A such that y(x) n D = 0. 
Indeed, for that to happen one would need that x = (x'-^'^^x'^-') with either x*^"*^' > 
(9/10)L and x^^) < -(9/10)L or x^^) < -(9/10)L and x^^) > (9/10)L, which is clearly 
incompatible with x G A, cf. (8.2). This shows that F(x^^^) is a simple curve. Observe 
that it is for this issue that it was important to change the shape of the domain from A 
to A. Of course, the value (9/10) in the definition of D could be changed to any number 
larger than 1/2. 
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Figure 5. The domain A, the sub-domain D (enclosed by the dashed hne; pro- 
portions are not respected) and the contour T{a) (thick hne) of a configuration 
(7 £ Q. Big full dots denote "vertex" sites (remark that they occur in nearest- 
neighboring pairs), empty dots denote "mountains" and * denote "valleys". Ver- 
tices can occur only in the right-most or left-most column which intersects M{a), 
or in the highest or lowest row which intersects M{a), so there are plainly at most 
8 vertices. The curve T{a) can be seen as the union of four monotone curves, 
which in the figure are delimited by the dotted lines. For each monotone curve, 
one has m(a) + w{a) — v{a) = 1. 



In all cases, the configuration belongs to G after the move, and the proof of a (t) E ^ for all 
t > is complete. 

To prove (8.7), assume that the last update before r^ consisted in hipping from "— " to "-I-" 
a "mountain" site x £ dD (if instead the move consisted in flipping from "— " to "-I-" a "vertex" 
site, the proof which follows would be very similar). Call x the spin configuration just before 
the last update and, just to fix ideas, assume that the two "— " neighbors of x in x are x — (0, 1) 
and X — (1,0). Since a; is a mountain and x £ ^j it is immediate to realize that the set 



Jix):=An{x + {i,j),i,j>0,i + j>0} 



(8.9) 



,W 



is a subset of A \ M{x), otherwise the condition |9r(x)| = '^'^i=i{urnax{x) ~ "^minix)) would 
be violated. Next, from the definition of the set D one sees that J(x) has cardinality at least 
ipL'^, uniformly in x G dD (explicitly, one can take ip to be slightly less than 2(1 ~ 9/10)^). The 
estimate (8.7) is then proven. 
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Figure 6. In the top drawing, an update consisting in the flip of the "vertex" 
spin Gx {x is marked by a dot) and of its neighboring vertex site from the value 

"— " to "+". In this case, UmLx decreases by 1 (while Umax and u^^^, i = 1, 2 are 
unchanged) and T{a) remains a simple curve, but its length decreases by 2. In 
the middle drawing, the effect of flipping from "— " to "+" a "mountain" ax- r(o') 

remains a simple curve and neither its length nor the values {umax,Umin\ 1=1,2 
vary. Similarly, the third drawing shows the effect of flipping a "valley" spin. In 
the three cases, only a portion of dT[(T) is drawn. 

Theorem 7 

D 
The simplification of considering a dynamics which evolves in the set Q is that for such 
configurations the numbers of valleys, mountains and vertices satisfy simple relations: 

Lemma 7. Given a G Q, letv{a) (resp. m{(T),w(a)) be the number of valleys (resp. mountains, 
vertices) in A. Then, m{a) + w{a) — v{a) = 4 and w{a) < 8. 

The proof of Lemma 7 is best explained through a picture, and therefore we refer to the 
caption of Figure 5. 

Now we can finish the proof of (8.1). Thanks to Theorem 7, we see that (8.6) follows if we 
have 

P(|M(a-(0))| - \M{a-{coL^))\ > i^L^) < e'^^ (8.10) 

By the exponential Tchebyshev inequality one has for A > 



P(|M(a-(0))| - \M{a~{coL'))\ > ijjL') < e"^^^ E 



gA(|M(a-(0))|-|M(a-(coL2))|) 



.11) 



-AVi-L^ 



e"^^^ 0(coL^) (8.12) 
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where </>(*):= B[f{a-{t))] and 

J(^) .^gA(|M(a-(0))|-|M(a)|)_ 

Note that one has <p{0) = 1 and 

d(t){t) 



where we recah that fi^ is the law of ^"(t) and C is its generator. By the very definition of 
the modified dynamics, if M{a) ^ (-D U dD) then C{a, a') vanishes for every a'. If instead 
M(a) 5 (D U dD), we know from Theorem 7 that we need only to consider the case a G Q. 
Therefore, for M{a) 5 (D U dD) one has (recall the discussion after Remark 6) 



>C/(ct) = fia) 



^('^) 7, + "1(0") + W[(7 



2 y ' ' V 2 y ' ' V 2 



.14) 



Now we choose A = 1/L. Since clearly there exists a constant c such that v{a),in{a),w{(T) < cL 
for every o" G ^, we have 



Cfia) < fia) 



-{m{a) + 2w{a)-v{a)) + ^ 



c' 



< jfi^), (8.15) 



where we used Lemma 7. Plugging this inequality into (8.13), we find 

^ < ^0«) (ai6) 

which implies (J){cqL'^) < e'^ '^'^ and, together with (8.11), 

P {\M{a'{0))\ - \M{a-{coL^))\ > i,L^) < e-^^^-^'^^'l (8.17) 

Choosing cq < ip/c' we obtain (8.10) and therefore (8.1). 



Eq. (8.1) 

□ 



9. Proof of Theorem 3 



9.1. Q{l/L) louver bound on the gap: a perturbative argument. Here we prove the lower 
bound gap > c/L for /3 > ClogL with C large enough. The result is particularly interesting 
in d = 2, in view of the matching upper bound in Theorem 3, but as we mentioned the proof 
works also in d = 3. Actually, we give the proof in the three-dimensional case, which is slightly 
more complicated. 

We introduce the matrix U := {f^(o', cj')}{o-,o-'enA}' unitarily equivalent to the matrix C (the 
generator (2.7)), as 

U{a,a') := yG^)C{a,a')^l=. (9.1) 

The spectrum of U coincides with the spectrum of C, so that we have to prove that the smallest 
non-zero eigenvalue of — C/ is lower bounded by c/L. Note that U is symmetric thanks to the 
reversibility condition TT{a)C{a,a') = ■K{a')C{a' ,a). 

Remark 7. The transformation (9.1) is the analogue of the inverse of the "ground-state trans- 
formation" which maps a Schrodinger operator of the form H = —A + V (which acts on L^(M'^)) 
into the operator 

-C = A-^HA = -A - — VV'o • V, 
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where ipo{-) is the ground state eigenfunction (we assume for definiteness that HipQ = 0, i.e. 
the ground state energy is zero) and the unitary operator A acts as {Af){x) = ipo{x)f{x) for 
/ G L^(M ). Note that C is the generator of a diffusion with drift. 

Looking at the definition of U and C, one immediately reahzes that 

U{a,a) = C{a,a) = -Y,C{a,a'). (9.2) 

• if cr' = a^^' for some x G A and H^{a) = H^{a'), then 

U{a,a') = C{a,a') = -; 

• if cr' = a^^^ for some x G A and \H^{a') - H^{a)\ > (hence > 4) then 



exp(-(/3/2)(g+(a)-g+(a0)) ^ ^.^z^. 
l + exp(-/3(i:f+(a)-i:f+(^0)) 



t/(cx, a') = ^-Pl-l/^/-il-A " --A '^ < ,-2/.. (9.4) 



• if a 7^ o"' and there exists no x such that a' = a^^\ then £((T, a') = C/(cr, a') = 0. 

If we write U = Uoo + R with Uoo = lim^^oo U we see that ah the matrix elements of R are 
smaller than 1/L^ if /3 satisfies (2.13) with C large enough. Since each row of R has at most 
|A| = 0{L^) non-zero elements, one sees easily that the spectral radius of R is 0{L^'^). 

We are therefore left with the task of proving that the smallest non-zero eigenvalue of —Uoo 
is larger than c/L. Thanks to formulas (9.2)-(9.4), we have that 

Uoo{a,a) = -i|{x G A : i7+(a(-)) = F+(a)}| - |{x G A : i7+(^(-)) < i7+(^)}| (9.5) 

and, for a ^ a' ^ 

Uooia,a') = \l '\ . ^' = "^^^ ^"^ HUa') = HU<ry, ^,,^ 

[ U otherwise. 

If we decompose ^\ into a finite number M of equivalence classes Cj, 1 < i < M, where two 
configurations belong to the same class iff they can be connected via a finite number of single 
spin-flips which do not change the energy HJ^{-), then f7oo(o', o"') = whenever o" G Cj,(t' G Cj 

with i ^ j. In other words, one can write Uoo in a block matrix form Uoo = (Bff^iUoo where 

C/oo = {Uoo{o',o'')}{a,a'eCi} ^^'^ ill particular, if S{Uoo) denotes the spectrum of Uoo, one has 

element is the all "+" configuration, one has S{Uoo ) = {0} (cf. (9.5) and observe that any spin 

fiip increases the en( 
for some positive c. 
Let us fix i > 1, 
eigenfunction. We want to show that 



^{Uoo) = U^^5(f/(X) ). It is clear that, if Ci = {+} denotes the equivalence class whose unique 

ser 
fiip increases the energy if o" = +). We need to prove that, for every i > 1, S{Uc^) C (— oo, —c/L) 

(i) 

Let us fix i > 1, let A = Aj be the smallest eigenvalue of —Uoo and 5 : Cj i— )■ M an associated 



A > I (9.7) 



with c independent of i. The key point is the following: 
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Lemma 8. For every r],r]' £ Ci and t > 0, one has 

e*^- (r/, r,') = P" (a^Ct) = ??'; r > t) (9.8) 

where P^ denotes the law of the Ising evolution {a^'{t)}t>o in A at (3 = +oo with + boundary 
condition, started from the configuration rj, and r is the random time 

T = inf{t > : H+{ait)) < Hj^ir^)}. (9.9) 

Lemma 8 

Proof. Just check that the time derivatives of left- and right-hand side are the same. D 

(i) 

Remark 8. Note that Uso is just the generator of the "killed" Markov process which coincides 
with the Ising evolution except that it is killed when it exits the set Ci, cf. (9.8). 

Since the matrix {exp(tC/o^) (7^,77')}^^^/ g^, has strictly positive entries for t > (this follows 
from (9.8) and the definition of the equivalence classes Ci), the Perron- Frobenius theorem implies 
that the eigenvalue A is non-degenerate and the eigenfunction g can be chosen strictly positive 
on Ci. Normalizing g so that ^(^g(j. 5(c)^ = 1, we have 

e"'*= Yl 5(^)e*^^'(a,a')5(^')- (9-10) 

The desired estimate (9.7) then follows from Lemma 8 by letting t — )■ 00, once we prove the 
following result: 

Theorem 8. There exists a{L) < 00 and c > independent of L such that, for every i > 1, 
rj, f]' G Ci and t > one has 

pv (^^v^t) = ri'-T>t)< a{L) exp (-c j\ . (9.11) 

Proof. Let V be the smallest parallelepiped which contains the set {x G A : r/^; = — }, and let 
V~^ be the rectangular layer of points of V with maximal vertical coordinate. Plainly, t < Ty 
where 

TV = inf{t > : al{t) = + for every x G F+j : (9.12) 

this is because when the last "— " spin in V^ flips to "-I-", the energy decreases by at least 8. 
Then, 

P^? (a'?(t) = ri'-T>t)< P'' (ry > t) < P*?^ {tv > t) (9.13) 

where rjv is the configuration where spins take value "— " in V and "-1-" in A \ y, and we 
used monotonicity of the dynamics (the event Ty > t \s decreasing, and rjy < rj). Again by 
monotonicity, we can assume that V is the entire domain A, and that spins in y \ V^ are frozen 
at the value "— " during the entire evolution: both these operations make Ty stochastically 
larger. But, in this case, V~^ is just a (2L -|- 1) x (2L -|- 1) square, and the evolution in V^ 
coincides with the evolution of the two-dimensional Ising model with -|- boundary conditions 
and /3 = -|-oo (the "-I-" boundary conditions on the upper face of V~^ and the "— " boundary 
conditions on the lower face compensate exactly). We have then 

P''^ (Ty >t)< P('^=2) (r+ > t) (9.14) 

where p("=2) jg ^j^g ]^^.^ q£ ^j^g evolution of the two-dimensional Ising model in the (2L -|- 1) x 
(2L -|- 1) square with "-I-" boundary conditions started from the "— " configuration, and t+ was 
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defined in (2.17) as the first time when all spins in the square are "+". Finally, one has 



,(d=2) 



T+>t)< 



p(^=2) > 



L 



ltc/L^\ 



< a{L)e'^^L (9.15) 



gap < _^^,^^ ^ inf ^ — -y-^T — , (9.16) 



where in the first inequality we used monotonicity (if t+ has not been reached at time i L/^ /c, i = 
1, . . . , [tc/L^J — 1, we restart the dynamics from the all "— " configuration) and in the second 

Theorem 8 

one the result (3.4) of Theorem 2. This ends the proof of Theorem 8. D 

Remark 9. The asymptotic behavior exp(— ct/L) can be understood by comparison with the 
symmetric simple exclusion process (SSEP), as already noticed in [7]. In fact, by monotonicity, 
the tail of the law of t+ is bounded from above by the tail of the hitting time tq defined as 
follows. Consider the SSEP on [— L, L] starting with all the negative sites occupied and all the 
non-negative ones empty and define tq to be the first time at which the site "L" is occupied. 
In turn, using beautiful results by Liggett and by Arratia [1], the tail of tq is controlled from 
above by the tail of the same hitting time but for the process in which the L particles evolve 
as independent random walks without the exclusion constraint. Finally the latter has a tail 
exp(— t/L) simply because the tail for a single random walk is exp(— t/L^) and there are L of 
them. 

9.2. Upper bound on the gap. We will prove here the upper bound gap < l/(cL) for the 
dynamics in the square box A = {— L, . . . , L}^ . We use the variational characterization (2.6) of 
the spectral gap and we note that, if y is some subset of the configuration space VL\, one has 

7r(3^=) f:f\yc=0 Trip 

where the infimum is taken over functions which vanish on y^: just observe that 

Var^(/) = \j{f{a)-f{r))\{da)7:{dr) (9.17) 

> \ jUi^) - f{r))Md<yMdr){Uey^ + l.^yc) = T^iy^^Mf). (9.18) 

The bound gap < l/(cL) is therefore proven if one exhibits a set y such that 7r(3^'^) > 1/2 and 
a function / which vanishes on y^ and such that 

n{p) - L- ^^-'^^ 

For a G Oa, we let M{a) and T{a) be as in Definition 7 and we let Q C M^ be the square 
[— L — 1/2, L + 1/2]^. We consider the eight points of dQ which are at distance L/2 from a 
corner of Q and we call them Pi,i = 1, . . . , 8, with the convention that we order them clockwise, 
starting from the right-most one on the north side of Q. Let us assume for definiteness that L 
is even. 

Definition 9. We let y be the subset of configurations in Da such that T(a) is a simple curve 
of length 4(2L + 1) and such that pi £ dT{a), i = 1, . . . ,8. 

Note that 4(2L + 1) is just the length of dQ. Also, if o" G 3^ then the following properties are 
immediately verified, cf. Figure 7 : 

(1) the portion of dT{a) which connects p2i to p2i+i,i = 1,2,3 or ps to pi is a straight 
segment of length L + 1 included in dQ; 
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v(4) 



P7f 



P6/^ 



M 



P& 



P5 






♦_--'' 



P4 



v(l) 



(,y 



-7-^-® V'i 



VI 



v(2) 



Figure 7. The domain A with L = 6 and a configuration cr G 3^: only the 
boundary dV{o^ is drawn (thick hne). The four portions 7(*',i = 1,...,4 of 
5r(cr) cannot intersect and, if suitably rotated, are just lattice paths. According 



to formula (9.21), one has in this example V_3 = V 



7(1) 



V 



(1) 



-1 and V 



(1) 



V 



(1) 



VI 



(1) 



+1. 



(2) the portion of dV(a^ which connects P2j-i to P2j,^ = l, • • • ,4 (call it 7(*'') is a "lattice 
path" of length L (in the sense of Section 6). More precisely, if 7'*' is rotated counter- 
clockwise by an angle 7r/4 + 7r/2(i — 1), then expanded by a factor v2 and finally is 
suitably translated, then it becomes the graph of a function x G [— L/2,L/2] 1— )■ 0^. := 

\(j) G [-L/2,L/2] such that 



#' 










V» := 






(0 



G {-1, +1} for X G Z n [-L/2, L/2 - 1] 



(\y is linear in each interval (x, x + 1), x G Z n [— L/2, L/2 — 1] 
The test function / in (9.19) is then 



(9.20) 
(9.21) 
(9.22) 



where 



A4>'^' 



/(^):=l{<xe3^}n^(<^^*^(^))' 



n 



i=l 



W (L/2)-l 
/7rx\ (,1-Va; )/2 
COS I 



-L/2 



V L 



n 

x=0 



COS 



vr(x + 1^ 

z 



(l+vi'')/2 



(9.23) 



(9.24) 
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with the convention that 0^ = 1. Note that / vanishes on y^ and that 

T^{y^) > vr(o-^ = + for every x G A) = 1 + o(l) > 1/2 

as required above. 

Let p be the uniform measure over the set of lattice paths (p satisfying (9.20)-(9.22), and let 
£, be the generator 

L/2-1 

^9{4>)= Yl [pi9\<Py,yez,y^x)-g{^)]- (9.25) 

x=-L/2+l 

Remark 10. Conditionally on o" G 3^, the four lattice paths 0'*^ are independent and, since all 
their allowed configurations have the same length, '7r(-|3^) is just the product uniform measure 
/j®^, where each copy of p acts on a path (f)^'^', i = 1, . . . , 4. Indeed, the Ising Hamiltonian (2.2) 
equals 

H+{a) = const + 2|r(cr)| (9.26) 

which does not depend on o" if o" € 3^, cf. Definition 9. 

Lemma 9. /// is defined as in (9.23), one has 

Proof. Let us consider first the denominator of the left-hand side: 

4 \ 



7r(/2)=vr(3;)vrm5' 
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y] =^{y)W)Y (9.28) 



where we used the fact that 7r(-|3^) = p®\ see Remark 10. 

The numerator requires more work. From the definition (2.7) of the generator one has 

4 

n{f{-C)f) = Y, vr(a) n5('A«(^)) E (/(^) " ""-^M)) (9-29) 

where tTz,(t{') is defined in (2.4). The first observation is that, if /? > ClogL with C large enough, 
then 

|vr.,.(/)-/(a)|<-l/(^) (9.30) 

whenever a ^ y and z is such that the configuration o"'^-* satisfies |r(o"''^-')| > |r(cj)| + 1 (in 
particular, o"'^' does not belong to 3^). The reason is that, since /(o"'^-*) = 0, one has 

T^zAf) - /(^) = -^zA-^z)f{a), (9.31) 

where i^z^ai—f^z) is the probability, under Tr^^o-, that the spin at z equals — Uz, and this probability 
is smaller than exp(-2/3) < L"^ by the assumption that |r(cr(^))| > |r(o-)| + 1 (cf. (9.26)). 
CaU B{a) the set of z G A such that |r(o-('^))| = |r(cr)|. One has then 

4 

i^{f{-c)f)<Y<'y)\{9{<P^'^{'y)) E U\<y)-^zAf)) + ^^<f) (9-32) 

aey i=l zeB{cr) 
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since, if a € y, there are no spin flips which decrease |r((T)| (i.e. the energy). Finally, recalling 
Remark 10 it is an easy task to check the identity 

4 

j;^(a)nff(0^^^(^)) E ifi<r)-7T.Af))='^<y)Wfpi9i-C)g), (9.33) 

aey i=l z£B{(t) 

the reason being that the flips which leave |r(cj)| unchanged just correspond to the local updates 
in the generator C of (9.25). The factor 4 is due to the symmetry among the lattice paths 

Lemma 9 

(/>«,z = l,...,4. D 

The proof of (9.19) is concluded once we prove 

Theorem 9. There exists c > such that 

pMz^ < f (9.34) 

p(r) L 

Proof. Using the form (9.25) of the generator and the fact that p is the uniform measure over 
lattice paths, one sees that 

^ 1 ,^ [cos(7ra;/L) - cos(7r(x - l)/L)]^ , 2-, n (r,o^\ 

P{9{-C)g) = - E cos^(vrx/L) P (<7 l{v._.=+i,v.=-i}) (9-35) 

x=-L/2+l ^ ' ' 

— 1 9 

1 ^-^ [cos(7r(a; - 1)/L) - cos(7ra;/L)] , ^ . 

+ 2^_Z.^_^^ cos2(4x-l)/L) P^^ l{V.-i=-i,V.=+i}) 

where we recall that V^; = 0x+i — ^x and we used the symmetry x -H- —x to restrict the sum to 
re < 0. Clearly, one has 

E [cos(7rx/L) - cos(7r(x - 1)/L)]^ < cjL. (9.36) 

X 

Moreover, by an "equivalence of ensembles" argument, one has 

^ (^'^iv.-.=+i.v.=-i}) ^ , ,^,2^^^,L^ (9.37) 

pir) 

and a similar estimate for p ((7^1|y^_^^_i y^=_|_i|) I p{g^) where c' is finite uniformly in x < —1. 
The statement of the Theorem then immediately follows from Eqs. (9.35), (9.36) and (9.37). 

The proof of (9.37) is a direct adaptation of the proof of [3, Prop. 3.8]. Loosely speaking, 
the idea is to observe that the law of \y x)x=-Lj2,...,Ll2-\ under p{- g"^) / pig^) is a product 
law of independent (but not identically distributed) random variables conditioned to the event 

J /2 1 

^ __Ti2^x = 0. Next, one shows that, when computing the average of a local function like 
l{Vj;_i=+i,Vi=-i}) the conditioning can be eliminated, at the price of losing a multiplicative 
constant. Under the unconditioned product law, the expectation of l{v^_i=+i,Vi=-i} is easily 
computed and turns out to be proportional to the r.h.s. of (9.37). 

Theorem 9 
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